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ABSTRACT 

Context. Asymptotic giant branch (AGB) stars are one of the largest distributors of dust into the interstellar medium. However, the 
wind formation mechanism and dust condensation sequence leading to the observed high mass-loss rates have not yet been constrained 
well observationally, in particular for oxygen-rich AGB stars. 

Aims. The immediate objective in this work is to identify molecules and dust species which are present in the layers above the 
photosphere, and which have emission and absorption features in the mid-infrared (IR), causing the diameter to vary across the N- 
band, and are potentially relevant for the wind formation. 

Methods. Mid-IR (8-13 jim) interferometric data of four oxygen-rich AGB stars (R Aql, R Aqr, R Hya, and W Hya) and one 
carbon-rich AGB star (V Hya) were obtained with MIDIA^LTI between April 2007 and September 2009. The spectrally dispersed 
visibility data are analyzed by fitting a circular fully limb-darkened disk (FDD). 

Results. The FDD diameter as function of wavelength is similar for all oxygen-rich stars. The apparent size is almost constant between 
8 and 10 fim and gradually increases at wavelengths longer than 10 fim. The apparent FDD diameter in the carbon-rich star V Hya 
essentially decreases from 8 to 12 ywm. The FDD diameters are about 2.2 times larger than the photospheric diameters estimated from 
K-band observations found in the literature. The silicate dust shells of R Aql, R Hya and W Hya are located fairly far away from the 
star, while the silicate dust shell of R Aqr and the amorphous carbon (AMC) and SiC dust shell of V Hya are found to be closer to the 
star at around 8 photospheric radii. Phase-to-phase variations of the diameters of the oxygen-rich stars could be measured and are on 
the order of 15% but with large uncertainties. 

Conclusions. From a comparison of the diameter trend with the trends in RR Sco and S Ori it can be concluded that in oxygen-rich 
stars the overall larger diameter originates from a warm molecular layer of H2O, and the gradual increase longward of 10 jum can 
be most likely attributed to the contribution of a close AI2O3 dust shell. The chromatic trend of the Gaussian FWHM in V Hya can 
be explained with the presence of AMC and SiC dust. The observations suggest that the formation of amorphous AI2O3 in oxygen- 
rich stars occurs mainly around or after visual minimum. However, no firm conclusions can be drawn concerning the mass-loss 
mechanism. Future modeling with hydrostatic and self-consistent dynamical stellar atmospheric models will be required for a more 
certain understanding. 
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1. Introduction 

Asymptotic giant branch (AGB) stars are among the most im- 
portant distributors of dust into the interstellar medium due to 
their high mass-loss rates in combination with an eflPective dust 
condensation. Especially, the dust plays a crucial role for the 



* Based on observations made with the Very Large Telescope 
Interferometer (VLTI) at the Paranal Observatory under program IDs 
079.D-0140, 080.D-0005, 081.D-0198, 082.D-0641 and 083.D-0294. 

Color versions of the figures and Tables A to C.2 are only avail- 
able in electronic form via http://www.aanda.org. FITS files of the cal- 
ibrated visibilities are only available at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz- 
bin/qcat?J/A-FA/vol/page. 

pgiiQ^ Qf International Max Planck Research School (IMPRS) 



formation and acceleration of the dense wind. By providing the 
seed particles for interstellar grains, AGB stars contribute to the 
chemical evolution of the interstellar medium (ISM) and facil- 
itate further star and planet formation. Progress in the theoreti- 
cal understanding has been made (Hofner et al., 2003; Woitke, 
2006a; Hofner, 2008; Norris et al., 2012), but the wind forma- 
tion mechanism and dust condensation sequence in oxygen-rich 
AGB stars need further observational constraints. Also for low 
mass-loss carbon-rich objects the wind formation has not been 
fully understood (Mattsson & Hofner, 201 1; Sacuto et al., 201 1). 

The slow wind of AGB stars is driven by stellar pulsation 
in combination with radiation pressure on dust. In order to sup- 
port this current understanding, the location and composition of 
newly formed dust as function of the pulsation cycle, mass-loss 
rate and underlying chemistry are investigated. The size of the 
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inner region free of dust is one of the parameters best constrained 
from mid-IR interferometric observations. This quantity is fun- 
damental for understanding the condition of dust formation. 

The low surface gravity and the stellar pulsation lead to an 
increased scale height of the atmosphere and make it possible 
to provide the conditions for dust grain formation. The conse- 
quence of a very extended atmosphere is a not clearly definable 
and measurable photospheric radius Rphot^ (e g- Mennesson et 
al., 2002; Tej et al., 2003; Ohnaka, 2004; Woodruff' et al., 2004; 
Fedele et al., 2005) and observations at diff'erent wavelengths 
probe diff'erent atmospheric layers (Baschek et al., 1991; Scholz, 
2001). The measured apparent diameters are correlated to the ab- 
sorption and emission features of the most abundant and radia- 
tively important molecular species (e.g. Hofmann et al., 1998; 
Jacob et al., 2000, and references therein) as well as first dust 
grain species with high sublimation temperatures (e.g. Lorenz- 
Martins & Pompeia, 2000; Verhoelst et al., 2009). 

The immediate objective in this work is to identify molecules 
and dust which are present in close layers above the photo- 
sphere (at around 2 Rphot)^ and which have absorption features 
in the mid-IR, causing the diameter to vary across the N-band. 
Candidates in 0-rich environments are H2O, SiO, CO, and TiO 
molecules, and dust grains composed of Mg2Si04, MgSi03, 
Si02, AI2O3, and Ti02 (e.g. Woitke, 2006a). In C-rich stars, 
molecular layers of C2H2 and HCN, and dust composed of SiC 
and amorphous carbon (AMC) can be expected. 

The size of the photosphere, close molecular and dust lay- 
ers are constrained by a broad sampling of the visibility in the 
first lobe, i.e. by carrying out observations over a large projected 
baseline range. Visibility points at the shortest baselines are used 
to measure the contribution of surrounding dust in the more ex- 
tended circumstellar environment. Observations were conducted 
at several orientations to establish whether the sources are elon- 
gated or not. 

From infrared spectroscopy, it is already known that there is 
no simple relation between the spatial distribution of diff'erent 
molecules and pulsation phase (e.g. Woodruff' et al., 2009, and 
references therein). By monitoring the stars regularly over a few 
pulsation cycles the dynamic behavior is studied, in particular 
the variation in the distribution of the close warm layers of dust 
and molecules. 

One star in this study (W Hya) was already extensively de- 
scribed in Zhao-Geisler et al. (2011) (thereafter paper I) and is 
included here only for completeness. In addition, some content 
refers to paper I. Section 2 gives an overview of the target stars, 
and Sect. 3 describes the observation and data reduction. This is 
followed by Sect. 4 showing the light curves, spectra and results 
of the visibility modeling. The data are interpreted and discussed 
in Sect. 5 including an investigation of the dynamic properties 
and morphology. A summary is given in Sect. 6. 

2. Target properties 

The observations concentrated on a relative small number of 
stars to obtain a fair number of visibility points to sample the uv- 
plane as well as the pulsation phase. In total five stars, namely 
R Aql, R Aqr, R Hya, W Hya, and V Hya, were observed in 
the framework of a Guaranteed Time Observation (GTO) pro- 
gram (cf. Sect. 3.1). The targets were chosen to cover different 
chemistries, evolution stages and mass-loss rates. R Aql, R Aqr, 
R Hya, and W Hya are oxygen-rich AGB stars, with R Aqr being 

^ The definition of /?phot = ^phot/2 used in this work is given in 
Sect. 4.4. 



a symbiotic system, while V Hya is an evolved carbon-rich AGB 
star with a fast coUimated wind. 

In the following, important theoretical and observational 
work from previous publications are summarized for each star 
giving an overview of their characteristic properties. Table 1 lists 
some of the relevant phenomenological features for each star. 
The given K band diameter, but also other optical and IR in- 
terferometric angular diameter measurements, are discussed in 
Sect. 4.4. The distances given in Table 1 are assumed through- 
out this paper. 

2.1. RAquilae 

The visual period of the Mira R Aql has declined quite dra- 
matically with an average rate of approximately 0.4 days per 
cycle since 1900 (Greaves, 1998; Greaves & Howarth, 2000). 
While in 1915 the period was about 320 days, it declined to 
about 264 days in 2010. This has been accompanied by a de- 
crease in amplitude of about 1 mag (Bedding et al., 2000). The 
shrinking of the period is attributed to a recent thermal pulse 
(Wood & Zarro, 1981). R Aql is listed in the Washington Visual 
Double Star Catalog (WDS Mason et al., 2001) as triple sys- 
tem. However, the two companion stars have dissimilar proper 
motions and are probably only optical companions (Greaves & 
Howarth, 2000). 

The distance to R Aql is relatively well known. The 
Hipparcos catalog (Ferryman & ESA, 1997b) gives a value of 
211^42 PC, while Kamohara et al. (2010) estimated the dis- 
tance to 214^32 pc via maser observations. The mass and lu- 
minosity, obtained from modeling and observations, are 1 Mq 
and 3470 L©, respectively (Hofmann et al., 2000). Pulsation 
phase dependent eff'ective temperature determinations range 
from (2550 ± 150) K (Haniff' et al., 1995), over 3072 ± 161 K 
(Hofmann et al., 2000) to 3198 ± 147 K (van Belle et al., 1996). 
From the relations between pulsation period, mass and linear ra- 
dius, Hofmann et al. (2000) came to the result that R Aql is a 
fundamental mode pulsator in agreement with recent considera- 
tions. 

SiO, H2O and OH maser emissions have been monitored for 
this 0-rich star. SiO maser emission was reported by Benson et 
al. (1990), Pardo et al. (2004) and Cotton et al. (2010). The lat- 
ter author derived for the 43.1 and 42.8 GHz emission average 
ring diameters of (26.1 ± 1.5) mas and (22.3 ± 1.9) mas, respec- 
tively, in agreement with the supposed extended atmosphere. At 
a larger diameter of around 330 mas. Lane et al. (1987) and 
Brand et al. (1994) detected H2O maser emission. Studies of 
OH masers were performed by Bowers et al. (1989), Etoka & 
Le Squeren (2000) and He et al. (2005), showing that this maser 
emission originates in a region with a diameter between 2000 
and 3000 mas. 

Matsuura et al. (2002) modeled the extended molecular 
spheres of R Aql with two layers of water vapor of diff'erent tem- 
peratures. The radius of the hot layer changes from 1 to 2 pho- 
tospheric radii between visual minima and maxima. The model 
could explain the emission and absorption features seen in the 
near infrared (2.5 - 4.0 ^m). While Cotton et al. (2010) could 
not identify strong asymmetries in the SiO maser distribution. 
Lane et al. (1987) and Bowers et al. (1989) saw highly elon- 
gated and complex structures in the H2O (NW-SE oriented) and 
OH masers (NE-SW oriented), respectively. In addition, the in- 
terferometric measurements by Ragland et al. (2006) gave a sig- 
nificantly non-zero closure phase, indicating an asynmietry as 
well. 
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Target name: 


RAql 

HIP 93820 


RAqr 

HIP 117054 


RHya 

HIP 65835 


WHya 

HIP 67419 


VHya 

HIP 53085 


Associated calibrator^: 


77 Sgr, 6 Peg 


30 Psc 


2Cen 


2Cen 


a Hya 


Type: 


Mira 


Mira 


Mira 


Mira/SRa 


SRa/L 


Chemistry: 


O-rich 


O-rich 


O-rich 


O-rich 


C-rich 


Mass'' (Mo): 


1.0 


1.5 


2.0 


1.0 


1.0 


Distance^ (pc): 


914+45 


250 + 50 


130 + 25 


QO+30 
^^-18 


360 + 70 


12 jum flux^ (Jy): 


402 + 36 


1580 ± 60 


1590 ± 80 


4200 ± 210 


1110 ±60 


Spectral classification^: 


M7 III e var 


M7 III pev 


M7 III e g 


M7 III e g 


N; C7,5e 


Effbctive temperature'' (K): 


3000 + 300 


2800 + 300 


2700 + 300 


2500 + 300 


2650 + 300 


Luminosity'^ (L©): 


3470 + 500 


5000 + 500 


11600+ 1000 


5400 + 500 


7850 + 500 


K-band diameter^ (mas): 


11.5 ± 1.4 


16.6 ± 3.0 


24.5 ± 1.2 


42.8 ± 3.2 


14.5 ± 0.3 


Mass-loss rate§ (lO"'^ Mgyr"^): 


6-35 


0.1-9 


0.1-2 


0.8-5 


10 - 600 


Multiplicity'': 


triple system 


SB 


wide binary 


wide binary 


wide binary 


Evolutionary stage'': 


TP-AGB 


TP-AGB 


TP-AGB 


early-AGB 


post-AGB 


Jets/fast wind'': 


no 


yes 


no 


no 


yes 


Asymmetry'': 


maybe 


yes 


no 


yes 


yes 


Maser'': 


SiO, H2O, OH 


SiO, H2O, OH 


SiO, H2O, OH 


SiO, H2O, OH 


CO? 


Period'': 


decreasing 


stable 


decreasing 


stable 


two periods 


Others'': 


recent 
He flash 


symbiotic 

system 
(Mira+WD) 


detached 
shell, recent 
He flash 




common 
envelope, 
rapid rotator 



Notes. ^ For calibrator properties see Sect. 3.1 and Table 3. '' See text for references, errors are approximations derived from the values found in 
the literature. ^ For R Aqr, R Hya and V Hya derived from the period-luminosity relation (cf. Whitelock et al., 2008), and for R Aql and W Hya 
obtained from maser measurements (Kamohara et al., 2010; Vlemmings et al., 2003, respectively). ^ IRAS flux (Neugebauer et al., 1984); 
^ SIMBAD; ^ Values are averages from the literature (cf. Sect. 4.4); § Range of current mass-loss rates as given in Sect. 2. For W Hya see e.g. De 
Beck et al. (2010), MuUer et al. (2008) and Justtanont et al. (2005). 



Total mass-loss rates were estimated to 8 x 10 ^ Moyr ^ and 
5.6 X IQ-^ Moyr-i by Gehrz & Woolf (1971) and Hagen (1982), 
respectively, whereas Knapp & Morris (1985) found a much 
higher value of 3.5 x 10"^ Moyr"^ The dust mass-loss rate was 
derived by Hagen (1982) with a value of 6 x 10"^ Moyr"^ . Wind 
velocities are measured in a range from 7 to 10 km/s (Bowers et 
al., 1989, and references therein). 



2.2. RAquarii 

R Aqr is the closest known symbiotic binary at a distance of 
about 250 pc (derived from the period-luminosity relation by 
Whitelock et al., 2008). This D-type (dusty) symbiotic system 
consists of a 1 .0 - 2.0 M© Mira variable and a 0.6 - 1 .0 M© white 
dwarf (WD) (Hollis et al., 1997; Boboltz et al., 1997; Tatebe et 
al., 2006; Gromadzki & Mikolajewska, 2009) that accretes mat- 
ter through a disk (e.g. Hollis et al., 2000). The R Aqr binary 
system was for the first time resolved in the continuum radio 
emission at 7 mm by Hollis et al. (1997). Hollis et al. (1997) 
and Gromadzki & Mikolajewska (2009) derived orbital solu- 
tions with a period of about 44 years, while Mcintosh & Rustan 
(2007) estimated the orbit period to be 34.6 years. The fitted 
semi-major axis is in all cases on the order of 15 AU (60 mas), 
meaning that even during periastron passage, the Mira variable 
remains relatively far from filling the Roche lobe (Gromadzki & 
Mikolajewska, 2009). 



The binary system hosts a compact Hii region within a fila- 
mentary oval nebula of 30 arcsec size (e.g. Kafatos et al., 1986). 
It is surrounded by a large and expanding hour-glass shaped neb- 
ulosity with an extension of at least 120 arcsec (Wallerstein & 
Greenstein, 1980; Hollis et al., 1985; Henney & Dyson, 1992; 
Corradi et al., 1999). The accretion disk, formed around the com- 
pact component, gives rise to prominent jets, detected across all 
spectral domains (cf. e.g. Kellogg et al., 2007; Nichols & Slavin, 
2009; Hollis et al., 1991; Paresce & Hack, 1994; Hollis et al, 
1985). The synmietric jets, oriented along a NE-SW axis, ex- 
tend up to 2500 AU (10 arcsec) and have expansion velocities of 
90 to 200 km/s. 

The current mass-loss rate of this star is rather low, but un- 
certainties are very high due to diflficulties in obtaining reliable 
gas-to-dust ratios. Rate estimates range from 1.3 x 10"^ Moyr"^ 
(Henney & Dyson, 1992) and 6 x 10"^ Moyr"^ (Matthews & 
Reid, 2007) to 3 x 10"^ Moyr'^ (Mennesson et al., 2002) and 
8.9 X 10"^ Moyr~^ (Danchi et al., 1994). In most calculations 
an eflfective temperature of 2800 K (Burgarella et al., 1992; 
Matthews & Reid, 2007) and a luminosity of 5000 L© (Ragland 
et al., 2008) were assumed. 

The O-rich Mira variable R Aqr is one of only three among 
48 symbiotic Miras that exhibits H2O masers (Ivison et al., 1994, 
1998; Whitelock et al., 2003). A complementary study of SiO 
masers showed strong emission at 42.8 and 43.1 GHz. SiO maser 
ring diameters were obtained by Boboltz et al. (1997), Hollis et 
al. (2001) and Cotton et al. (2004, 2006), with relatively con- 
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sistent values between 30 and 33 mas and errors on the order 
of 1.5 mas. This corresponds to a radius of 1.9 photospheric 
radii. OH maser and CO thermal lines are very weak and only 
tentative detections have been published by Ivison et al. (1998) 
and Groenewegen et al. (1999), respectively. In the most pop- 
ular model, UV radiation and a fast wind from the companion 
remove the outer envelope of dusty molecular gas, where an OH 
maser or a thermal CO line could originate (Ivison et al., 1998). 

Ragland et al. (2008) modeled the star with a three- 
component model, consisting of a synmietric central star sur- 
rounded by a water vapor shell with a radius of about 2.25 pho- 
tospheric radii, and an off-axis compact feature at about 2 photo- 
spheric radii at a position angle of 131°. They concluded that the 
observations are best explained with a clumpy, extended H2O 
circumstellar envelope. In this model the SiO masers appear at 
the outer edge of the molecular envelope, as reported for other 
Mira stars. 

Possibly caused by interactions with the close companion, 
asymmetries could be identified in the close environment around 
R Aqr. Ragland et al. (2006, 2008) found non-zero closure 
phases, and HoUis et al. (2001) and Cotton et al. (2004, 2006) 
detected an asymmetric distribution of the SiO maser emission 
with position angles between 150° and 180°. Terminal wind ve- 
locity measurements are rare. Kotnik-Karuza et al. (2007) ob- 
tained a value of 6 - 7 km/s, while Dougherty et al. (1995) as- 
sumed 10 - 30 km/s. 

2.3. R Hydrae 

The 0-rich variable star R Hya is well known for its declining 
period (Wood & Zarro, 1981; Zijlstra et al., 2002) and the pres- 
ence of a detached shell. In the past few centuries, the period of 
R Hya has declined by over a hundred days and has remained 
constant since 1950 (Zijlstra et al., 2002). The steady decrease 
in the period can be possibly attributed to a recent thermal pulse 
(Wood & Zarro, 1981; Zijlstra et al, 2002). In the post thermal- 
pulse evolution the decline in luminosity causes a reduction in 
stellar radius, which in turn, causes the period to become shorter. 

The detached shell observed around R Hya (Young et al., 
1993; Hashimoto et al., 1998) indicates a change in the mass- 
loss rate some 220 years ago. The shell is located about 1.9 arc- 
sec from the star 250 AU). While the mass-loss rate be- 
fore 1770 is estimated to be between 1.5 x 10"^ Moyr"^ and 
3 X 10"^ Moyr-i (Hashimoto et al., 1998; Zijlstra et al., 2002), 
the present day mass-loss rate is a factor of ^ 20 lower and be- 
tween 1 X 10"^ Moyr"! and 4 x 10"^ Moyr"^ (Decin et al., 2008; 
Teyssier et al., 2006, respectively). De Beck et al. (2010) gives 
an average mass-loss rate of 1.6 x 10"^ Moyr"^ 

This derived mass-loss history nicely agrees with the period 
decline analyzed by Zijlstra et al. (2002). Even the stellar evolu- 
tion tracks calculated by Vassiliadis & Wood (1993) confirmed 
that mass-loss fluctuations during the thermal pulse cycle can 
lead to detached circumstellar shells. An apparent large detached 
shell may also arise from the interaction of the AGB wind with 
the ISM. Ueta et al. (2006) and Wareing et al. (2006) detected a 
far-IR nebula at a distance of about 100 arcsec to the west. This 
is explained by a slowing-down of the stellar wind by surround- 
ing matter. Therefore, no extra mass-loss modulation at around 
100 arcsec, i.e. around 10,000 years ago, needs to be invoked 
(Decin et al., 2008). 

Distance estimations for R Hya range from (110 ± 21) pc 
(Jura & Kleinmann, 1992) to 165 pc (Zijlstra et al., 2002), while 
the period-luminosity relation gives an intermediate value of 
130 pc (Whitelock et al., 2008). The mass and luminosity, in- 



ferred from modeling and observations, are 2 M© and 1 1,600 L©, 
respectively (Zijlstra et al., 2002). Eff'ective temperature de- 
terminations range from 2600 K (Teyssier et al., 2006), over 
2680 ± 70 K (Haniff-et al., 1995) to 2830 K (Zijlstra et al., 2002). 
R Hya is believed to be a wide binary system with a very long 
orbital period. The angular separation is 21 arcsec (WDS). 

The wind in the inner CSE starts with a velocity of 1.5 km/s 
and accelerates to 6.5 km/s farther out, as modeled by Teyssier 
et al. (2006), before reaching a terminal velocity of 7.5 km/s 
(Justtanont et al., 1998) to 10.0 km/s (Hashimoto et al., 1998; 
Knapp et al., 1998). Justtanont et al. (1998) found that a 13 fim 
dust feature and the appearance of strong emission lines of CO2 
originate in a warm layer close to the star. R Hya exhibits maser 
emission of SiO (Humphreys et al., 1997), H2O (Takaba et al., 
2001) and OH (Lewis et al., 1995). The resolution of these ob- 
servations were too low to derive any spatial information. Thus, 
asymmetries were not studied with masers. A morphology study 
in the visual (0.7 - 1.0 //m, Ireland et al., 2004a) and K-band 
(Monnier et al., 2004) gave no indications of a departure from 
symmetry. 

2.4. W Hydrae 

The main characteristics of the 0-rich star W Hya were exten- 
sively described in paper I, but are summarized in Table 1. 

2.5. V Hydrae 

V Hya is a classical (N type) carbon star (Zuckerman et al., 
1977) and is believed to be in a short-lived but critical stage 
in the evolution of a mass-losing AGB star into a bipolar PN 
(Tsuji et al., 1988; Kahane et al., 1996). This dust-enshrouded 
star has a C/0 ratio of 1.05 and an eff'ective temperature of about 
2650 K (Lambert et al., 1986). The luminosity is estimated to 
be 7850 L© (Knapp et al., 1997). Distance calculations range 
from 340 pc (Barnbaum et al., 1995) and 380 pc (Knapp et al., 
1997) to 440 pc (Olivier et al., 2001) and 550 pc (Bergeat et 
al., 1998), while the period-luminosity relation (Whitelock et al., 
2008) gives a value at the lower end with 360 pc. 

V Hya has two variability periods, (529.4 ± 30) d with 
a peak-to-peak variation of 1.5 mag and (6160 ± 400) d 
(^11 years) with a peak-to-peak variation^ of 3.5 mag, and is 
classified as semi-regular (SRa) by Mayall (1965) and Mira vari- 
able by Knapp et al. (1999). While the 529.4 d period is typical 
of a luminous AGB star, Knapp et al. (1999) suggests that the 
regular long period dimming of V Hya is due to a thick dust 
cloud orbiting the star. 

The rotation velocity of an AGB star, evolving in isolation, is 
not likely to be greater than 2 kms~^ even if its main-sequence 
progenitor rotated at breakup velocity. In contrast, the rotation 
velocity, v sin /, derived for V Hya from a high resolution optical 
spectral broadening analysis by Barnbaum et al. (1995) is on the 
order of 11 to 14 kms~^ Barnbaum et al. (1995) and Kahane et 
al. (1996) concluded that this rapid rotation is due to the spin- 
up by a companion in a conmion envelope configuration, but 
alternative explanations have also been proposed (Luttermoser & 
Brown, 1992; Olivier & Wood, 2003). The secondary star could 
be an early G or F star, implying a mass of 1.0 - 1.5 M©, or 
possibly a WD (Barnbaum et al., 1995). V Hya is in contrast 
inferred to have a mass of about 1 M© (Kahane et al., 1996). 



^ The long secondary pulsation amplitudes in JHKL are 2.4, 2.1, 1.7, 
and 0.7 mag, respectively (Olivier et al., 2001). 
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In addition, a wide companion, with a separation of 46 arcsec, 
exists (WDS). 

A consequence of the fast rotation is an enhanced equatorial 
mass loss producing a disk and a jet-like structure (e.g. Soker, 
1992; Morris, 1987). High angular resolution CO maps, mil- 
limeter, infrared and optical spectra suggest that the circumstel- 
lar structure of V Hya consists of three kinematic components: a 
low-velocity disk with a radial velocity of Av » ± (8 to 16) km/s, 
a intermediate- velocity wind with Av ± 60 km/s, and a high- 
velocity jet with Av ^ ± (60 to 200) km/s (Zuckerman & Dyck, 
1986; Sahai & Wannier, 1988; Lloyd Evans, 1991; Knapp et al., 
1997; Olivier et al., 2001; Hirano et al., 2004; Sahai et al., 2009). 

The low-velocity circumstellar environment component is 
flattened or has a disk-like shape, and is elongated along the 
north-south direction (Tsuji et al., 1988; Kahane et al., 1988, 
1996; Sahai et al., 2003). This may enable or enhance the forma- 
tion of an accretion disk and supports a model in which the jet is 
driven by an accretion disk around an unseen, compact compan- 
ion. Images obtained at 9.8 and 1 1.7 jum by Lagadec et al. (2005) 
show an additional slightly elongated structure in the south-west 
direction, tracing the dust emission from material blown away, 
while the overall structure is roughly spherically symmetric. 

The star has a large infrared excess and strong millime- 
ter molecular line emission, showing that it is losing mass 
at a fairly high rate. Total mass-loss rates were estimated to 
6.1 X 10"^ Moyr-i (De Beck et al., 2010), 2.5 x 10"^ Moyr'^ 
(Knapp et al, 1997), 1.5 x 10"^ Moyr'^ (Kahane et al., 1996; 
Knapp et al., 2000), and 1.0 x 10"^ Moyr'^ (Bambaum et al., 
1995) with dust mass-loss rates of 2.0 to 5.7 x 10"^ Moyr'^ 
(Knapp et al., 1997; Hirano et al., 2004, respectively). Knapp 
et al. (1997) modeled the dust envelope of V Hya by assum- 
ing grains consisting of amorphous carbon with dimensions of 
0.2 fim. Zuckerman & Dyck (1986) discovered the presence of 
a narrow CO emission feature superposed on a standard broad 
stellar CO profile, which probably represents the first example of 
a CO maser ever seen in any interstellar or circumstellar source. 



3. Observations and data reduction 

3.1. Interferometric observations with MIDI/VLTI 

The data presented here were obtained with the mid-IR (8 - 
13 yc/m) interferometer MIDI (Leinert et al., 2003, 2004) at the 
Very Large Telescope Interferometer (VLTI) in service mode us- 
ing the Auxiliary Telescopes (ATs). All five stars were moni- 
tored from P79 to P83 (April 2007 to September 2009) under 
the program IDs 079.D-0140, 080.D-0005, 081.D-0198, 082.D- 
0641 and 083.D-0294 in GTO time. An overview of the course of 
observations is given in Table 2. It should be noted that observa- 
tions were still ongoing at the time of analyzing the data, and that 
observations not executed in a specific semester are shifted to the 
next semester. A complete observation log is given in Table C.2 
in the appendix. Projected baselines range from 11 to 71 m and 
the position angles (PA, ?^; East of North) are differently dis- 
tributed for each star. The uv-coverages are shown later in the 
left hand panels of Fig. 6. 

Before or after each target observation a calibrator star is 
observed with the same setup in order to calibrate the visibili- 
ties and fluxes. The properties of the calibrator stars are listed in 
Table 3. Necessary input parameters for the calibration are the 
angular diameter (model diameter from Verhoelst, 2005)^ and 



^ http://www.ster.kuleuven.ac.be/~tijl/MIDI_calibration/mcc.txt; see 
also ESQ CalVin database: http://www.eso.org/observing/etc/ 



the IRAS^ 12 jim flux. The second row lists the associated tar- 
gets for which the calibrator is used. The angular separation be- 
tween calibrator and target is given in brackets and shows that 
they are located relatively far away from each other as a conse- 
quence of the low number of available calibrators. However, this 
is fortunately not a big concern in the mid-infrared. In particu- 
lar, the visibility calibration is not aff'ected by difl'erent airmasses 
since both, the correlated and uncorrected flux, changes propor- 
tionally with it. Nevertheless, observations need to be carried 
out under good seeing conditions (typical < 1.4'') and a clear 
sky. Both are necessary to have a long coherence time to obtain 
fringes and to reduce the infrared background. 

In order to identify molecules and dust species, the mid- 
infrared fringes are spectrally dispersed. The high spectral res- 
olution GRISM mode of MIDI, with R = A/AA = 230, would 
be therefore of advantage, but introduces additional problems in 
the reduction process (e.g. the photometric channels are unfa- 
vorably illuminated). Hence, mainly the PRISM with a spectral 
resolution of 30 was used. Since the targets are bright enough, 
observations were executed in SCI -PHOT mode, where the pho- 
tometric and the interferometric spectra are recorded simultane- 
ously. This has the advantage that the photometric spectrum and 
the fringe signal are observed under the same atmospheric con- 
ditions. 

3.2. MIDI SCI-PHOT data reduction 

The standard MIA+EWS^ (version 1.6) data reduction package 
with additional routines for processing SCI-PHOT data (Walter 
Jafl'e, private communication) was used. Measurements at wave- 
lengths beyond 12.0 jim were excluded due to too low fluxes of 
the calibrator stars in that wavelength regime and therefore diffi- 
culties to determine the signal in the presence of a high infrared 
background^. The remaining wavelength range from 8 to 12 //m 
is then binned into 25 wavelength bins. A detailed description of 
the reduction process and how the errors are derived are given in 
paper I and Zhao-Geisler (2010). In the end, 32 of 40, 26 of 39, 
64 of 72, 75 of 83 and 48 of 63 observations could be adequately 
reduced for R Aql, R Aqr, R Hya, W Hya and V Hya, respec- 
tively. Rejected are observations where the reduction process 
failed or unphysical visibilities arose due to bad environmental 
conditions. 



4. Light curves, spectra and visibility modeling 
results 

4.1. Light curves 

In order to assign a pulsation phase to the observations, vi- 
sual data from the American Association of Variable Star 
Observers (AAVSO)^ and the All Sky Automated Survey (AS AS 
Pojmanski et al., 2005)^ are used. After binning the AAVSO data 



http://irsa.ipac.caltech.edu/Missions/iras.html 
^ http : //www. strw.leidenuniv.nl/~nevec/MIDI 

^ Remaining flux, coming from short term variations of the high sky 
background and instrumental imperfections, not subtracted by the un- 
chopping has to be determined and removed. The reduction routines 
have to distinguish between sky, tunnel and source flux, and have to 
know the slit position. This changes with AT position, projected base- 
line and pointing. If the flux of the calibrator is too low, its flux can not 
be reliably recovered and the calibration process fails. 

^ http://www.aavso.org/ 

^ http://www.astrouw.edu.pl/asas/ 
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Table 2. Number of visibility measurements used in this study (see Appendix C for a detailed observation log and paper I in the case of W Hya). 





Semester 


Program ID 


Disperser 


RAql 


R Aqr 


RHya 


WHya 


VHya 


79 


Apr, 2007 - 30* Sep, 2007) 


079.D-0140 


GRISM 


11 


10 


18 


22 


8 


80 


Oct, 2007- 3 Mar, 2008) 


080.D-0005 


PRISM 


6 


10 


16 


17 


26 


81 


Apr, 2008 - 30* Sep, 2008) 


081.D-0198 


PRISM 


15 


7 


16 


18 


8 


82^ 


Oct, 2008- 3 Mar, 2009) 


082.D-0641 


PRISM 





4 


12 


14 


20 


83^ 


Apr, 2009 - 30* Sep, 2009) 


083.D-0294 


PRISM 


8 


8 


10 


12 


1 



Notes. ^ The observation program was still ongoing at the time of analyzing the data. 



Table 3. Calibrator properties. 



Calibrator name: 
Associated target^: 

Spectral classification'': 
Parallax^ (mas): 
12 yL/m flux'^ (Jy): 
Angular diameter^ (mas): 
Quality flag^: 



HD 167618 



6 Peg 

HD 206778 



30 Psc 

HD 224935 



2Cen 

HD 120323 



a Hya 
HD 81797 



RAql (46.4°) RAql (42.8°) R Aqr (10.3°) WHya (6.1°) V Hya (23.8°) 

RHya (12.0°) 



M2 III 
21.87 + 0.92 

214 + 9 
11.66 + 0.04 
1 



K2 lb var 
4.85 + 0.84 

104 + 5 
7.59 + 0.05 
1 



M3 III 
7.86 + 0.94 

87 + 4 
7.24 + 0.03 
1 



M5 III 
18.39 + 0.74 

255 + 26 
13.25 + 0.06 
2 



K3 III 
18.40 + 0.78 

158 + 10 
9.14 + 0.05 
1 



Notes. ^ For target properties see Table 1. The angular separation between calibrator and target is given in brackets. ^ Hipparcos (Perryman & 
ESA, 1997b); IRAS flux and model diameter from Verhoelst (2005), http://www.ster.kuleuven.ac.be/~tijl/MIDI_calibration/mcc.txt; ^ ESO 
CalVin database, http://www.eso.org/observing/etc/ 



Table 4. Fitted visual light and mid-IR light curves^. 



visual (AAVSO+ASAS) 
Target P To Vo Vampi 

days JD mag mag 



A0 



mid-IR (MIDI) 

Jy 



-^ampl 

Jy 



RAql 
R Aqr 
RHya 
WHya 
VHya'' 



271 +4 
390 + 5 
374 + 5 
388 + 5 
529 + 10 



2452910 
2452480 
2453231 
2452922 
2450045 



8.30 
8.75 
6.95 
7.56 
8.50 



1.80 + 0.15 
2.75 + 0.20 
1.95 + 0.15 
1.41+0.10 
1.50 + 0.10 



-0.21 + 0.50 
+0.15 + 0.05 
+0.08 + 0.05 
+0.15 + 0.05 
+0.29 + 0.05 



377 + 50 
1080 + 100 
1930 + 200 
4910 + 200 
1330 + 200 



80 + 50 (21%) 
330 + 50 (31%) 
240 + 50 (12%) 
510+100(10%) 
580 + 200 (44%) 



Notes. ^ V(t) = Vampi sin [27r(t - To)/P] + Vo and N(t) = Nampi sin [27r(t - To)/P - A0] + Nq, respectively. The phase shift A0 gives the positive 
off'set in respect to the visual phase. ^ V Hya has a second period with P2 = (6160 + 400) days, To^2 = 2446937 days, Vo,2 = 8.50 mag and 
Vampu = (3.50 ± 0.50) mag (Knapp et al., 1999). 



into 10 day bins a simple sinusoid is fitted to the AAVSO and 
ASAS data over a period of about 10 years (2000 - 2010). The 
resulting periods P, Julian Dates of maximum brightness Tq (de- 
fined as phase 0.0), mean visual magnitudes Vo and visual semi- 
amplitudes Vampi are given in Table 4. The errors given are esti- 
mates derived from the fitting process. The uncertainties for Tq 
and Vo are of the same order as the uncertainties for the period 
and the visual semi-amplitude, respectively. Note that V Hya has 
two periods and the longer period is taken from Knapp et al. 
(1999). 

For each star the photometric data as well as the fit are plot- 
ted in Fig. 1 versus time and versus visual phase, covering the 
period of the MIDI observations. In addition, the mid-IR fluxes 
in Jy (averaged between 11.5 and 12.5 jim) obtained with MIDI 
are shown in the right hand panels as function of visual phase. 
As described in paper I, the MIDI fluxes are strongly error-prone 
due to difficulties in the reduction removing instrumental im- 



perfections and the large sky background which varies on short 
time scales. Additional scatter in the N-band flux originates from 
cycle-to-cycle variations caused by the movement of the dust 
shell in the case of a dust driven wind (e.g. Nowotny et al., 2010). 

However, a clear phase dependence is detectable in each of 
the phase-folded plots. The results of a sinusoidal fit are sum- 
marized in Table 4. The table gives the relative phase shifts A<p 
of the mid-IR maximum with respect to the visual maximum, 
the mean N-band fluxes A^o and the semi-amplitudes A^ampi with 
the corresponding percentage flux variations given in brackets. 
R Aql shows only a very small variation in the mid-IR by eye, 
suggesting that the amplitude is overestimated by the fit, and it 
is even consistent with no variation at all. The poor fit could be 
the reason why a mid IR maximum before the visual maximum 
is obtained in comparison to the other objects. 

The Mira R Aqr and the semiregular variable V Hya show 
the highest mid-IR flux variations. Both systems are close bi- 
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Fig. 1. Left: Visual light curves covering the period of the MIDI observations. Simple sinusoidal fits are included to determine the pulsation phases 
used throughout this paper. The MIDI fluxes at around 12 yum are shown as well with the magnitude scale given on the right. Right: Same as left, 
but plotted versus visual phase. The MIDI fluxes at around 12 yt/m are shown in the lower panels of each plot. The data within the shaded regions 
in both columns are used for size and flux variation studies in Sect. 5.3. The corresponding plots for W Hya can be found in paper I (Fig. 3 and 4, 
respectively). 



nary systems and might therefore contain a large amount of cir- 
cumstellar material. The fit to the mid-IR flux data of R Aqr is 
fair and the flux variation of about 30% seems to be real. Even 
though the mass-loss rate is low, the orbiting dust in this sym- 
biotic system could be the reason for an increased mid-IR flux 
variation. 



In contrast, the mid-IR flux of V Hya varies not only with 
a high amplitude but also with a large scatter, resulting in a fit 
which is not very robust. Since V Hya displays two superim- 
posed visual periods, an assignment of a visual phase is not an 
adequate description. Even if the short period is Mira-like, the 
cause and influence of the long secondary period is not well 
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known. Many different causes probably contribute to its large 
and varying mid-IR flux. In C-rich stars the production of car- 
bonaceous dust is more eff'ective, and hence more dust can radi- 
ate. The fact that V Hya is probably a common envelope system 
might be related to this. However, V Hya's evolutionary status 
as post AGB star, exhibiting the superwind phase with a fast and 
dense mass loss, might be the main reason for the varying mid- 
IR flux. 

R Hya had a very smooth and regular visual light curve over 
the last decade. MIDI observations were only possible around 
minimum visual light, since the period is close to one year. As a 
result, the mid-IR light curve could not be well constrained and 
the observed flux variation is relatively small, only on the order 
of 10%. The mid-IR and visual light curves of W Hya are very 
stable with only moderate amplitudes as well. For both stars, the 
low mid-IR flux variations might be related to the low mass-loss 
rates. Since W Hya and R Hya are fairly nearby, both stars have 
high absolute mid-IR fluxes. 

Even if the mid-IR flux variations of all five AGB stars are 
much smaller than the visual ones, they are still on the order of 
10% to 30%. Similar flux amplitudes in the mid-IR have been 
reported for other AGB stars as well (cf. references given in 
paper I). Except for R Aql, the mid-IR maximum occurs al- 
ways after the visual maximum at an average visual phase of 
0.15 ± 0.05. This phase shift is consistent with previous stud- 
ies of AGB stars (cf. e.g. Lattanzio & Wood, 2004; Smith et al., 
2006; Nowotny et al., 2010, and references therein). The reason 
for this phase-lag is most probably related to the dynamic pro- 
cesses of shock front and dust formation. 

Since the mid-IR fits are relatively uncertain, all data are in- 
terpreted with respect to the visual light curve in the following. It 
should also be kept in mind that folding consecutive cycles into 
one cycle might not always be appropriate, since the pulsation is 
not strictly regular, in particular for V Hya (see Fig. 1). 

4.2. Spectra 

All calibrated MIDI spectra for each star are median averaged 
over all phases and cycles and are shown in the spectral en- 
ergy distribution (SED) plots in the left hand panels of Fig. 2. 
The uncertainties are given by the standard deviation. The flux 
levels of ISO^ and IRAS^^ have been adjusted to coincide with 
the 12 fim flux. Photometric data from 2MASS and IRAS are 
plotted as well. A blackbody curve is over-plotted as guidance, 
assuming a diameter and temperature which best represents the 
data points by eye. However, due to the infrared excess, strong 
metallic oxide lines, molecule absorption and dust extinction, it 
is not expected that a blackbody curve fits the spectral data in 
an appropriate way. Individual spectra of the MIDI observations 
are given in the right hand panels of Fig. 2, showing the same 
features as the averaged MIDI spectra in the left hand panels. 

The ISO spectra of the 0-rich stars are mainly dominated by 
absorption bands of H2O between 2.5-3.0 yum (stretching mode) 
and 5.0 - 8.0 fim (bending mode), and an SiO absorption band 
between 8 and 9 jim (y = 1 - 0). Distinct absorption lines of CO 
at around 2.4 ^m, OH at 2.9-4.0 fim, CO2 at 4.25 fim and SO2 at 
7.4 fim can be seen in the spectra of some of the stars. From tem- 
perature investigations by e.g. Justtanont et al. (2004) for W Hya 
it could be derived that these absorptions originate from diff'erent 
molecular layers located in the close atmospheric environment 
of the star. OH and CO absorption bands arise mainly from a 

^ Sloan et al. (2003) and http://iso.esac.esa.int/ida/ 
http://irsa.ipac.caltech.edu/Missions/iras.html 



hot (about 3000 K), dense region very close to the stellar pho- 
tosphere, where H2O is still photodissociated by shocks. H2O 
and CO2 absorption bands originate from a layer with a temper- 
ature of 1000 K, i.e. a molecular layer (MOLsphere) farther out. 
The SiO molecule absorption arises in the same region where 
the H2O molecular shell exists and where SiO is still not bound 
in dust grains. 

All 0-rich stars in the sample show dust emission in their 
ISO spectrum with diff'erent peculiarities. The features between 
10 and 20 fim are a combination of emission from amorphous 
silicates at 9.7 ^um, amorphous AI2O3 at around 11.5 ^um, and 
MgFeO at around 19 //m. Regarding the 13 //m emission fea- 
ture, seen in R Aql, R Hya and W Hya, it is still under debate 
whether it comes from spinel (MgAl204) (e.g. Posch et al., 1999; 
Fabian et al., 2001; Heras & Hony, 2005) or corundum (crys- 
talline AI2O3, a-Al203) (e.g. Onaka et al., 1989; DePew et al., 
2006). W Hya is the only star which clearly shows the molecu- 
lar SO2 emission feature in the spectrum (Justtanont et al., 2004; 
Zhao-Geisler et al., 2011), while R Aqr is the only star which 
is completely dominated by amorphous silicate dust emission. 
In addition, R Aqr has an overabundance of SiO molecules as 
discussed in Angeloni et al. (2007). All stars also exhibit a re- 
markable infrared excess. In the dust emission scheme of Sloan 
& Price (1998b) R Hya is classified as SE2t (broad oxygen-rich 
dust emission), R Aql as SE5 (structured silicate emission), and 
R Aqr and W Hya as SE7 and SE8, respectively (classic narrow 
silicate emissions). 

If the individual and averaged MIDI spectra are compared 
with the ISO spectra (cf. Fig. 2), it becomes obvious that the sili- 
cate emission at around 9.7 jim is not detected except for R Aqr. 
This behavior can be attributed to instrumental characteristics. 
ISO has a much larger field of view (FoV) compared to MIDI. 
With a small FoV of about 1 to 2 arcsec^^ the emission of the 
extended silicate dust shell could therefore not be observed with 
MIDI. 

This applies to the nearby stars R Hya and W Hya. R Aql 
and R Aqr are located farther away and have nearly the same 
distance. Therefore, both stars should both either show or not 
show the silicate dust emission feature, but this apparent contra- 
diction can be solved. In Sect. 5.1 it will be seen that the dust 
shell around R Aqr is much closer due to the fact that R Aqr is a 
symbiotic system probably containing a large amount of dust. 

However, this non-detection allows to derive a lower limit 
for the inner boundary of the silicate dust shell. Assuming a 
conservative value of the FoV of 1 arcsec, most of the silicate 
dust emission originates from a shell with an inner diameter 
larger than about 100 (>220 AU), 49 ^phot (>130 AU), 

and 28 ^phot (>90 AU) for R Aql, R Hya, and W Hya, respec- 
tively. However, this does not mean that silicate dust does not 
exist closer to the star. Its abundance is just below the detection 
limit of MIDI. For R Aqr it can be speculated that silicate dust 
must exist farther out, at a region with a diameter larger than 
72 

^phot (>250 AU), since with the larger FoV of ISO more sil- 
icate emission could be revealed than with the smaller FoV of 
MIDI (Fig. 2c). 

The absolute flux levels for the MIDI spectra are lower than 
those of ISO/SWS and IRAS as it should be expected since 
the FoV of MIDI is smaller than those of ISO/SWS and IRAS. 
Also the 3.5 fim region of the ISO/SWS spectra, considered as a 
pseudo continuum, does nicely overlap with the blackbody con- 



The exact value depends on the baseline used, the slit/mask position 
and other instrumental specifications. 

Photospheric diameter of the star, cf. Sect. 4.4. 
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Fig. 2. Left: Spectral energy distribution for all stars (see text). The black body curve is only included for guidance. Right: All available individual 
MIDI spectra around the visual maxima for all stars (including an average of these) except for R Hya were the spectra around the visual minima 
are given. The visual light phases are given in brackets. Note that due to the telluric ozone feature the uncertainties at around 9.6 yt/m are increased. 
The corresponding plots for W Hya can be found in paper I (Fig. 5). 



tinuum. However, this is not the case for R Hya. This might be 
due to a calibration problem at around 9 yum since the IRAS spec- 
trum is flat around that point while it is not for the ISO/SWS one. 
At the same time, the flux level of the MIDI spectrum is slightly 
higher, but this is probably due to the large MIDI flux calibra- 



tion uncertainties. A flux variation due to the pulsation can be 
ruled out since the ISO spectrum is taken at around visual maxi- 
mum while the averaged MIDI spectrum of R Hya is mainly an 
averages of spectra taken around visual minimum. 
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V Hya is the only carbon star in the sample. The MIDI spec- 
trum shows very clearly the dust emission of SiC at 11.3 jjm. 
Amorphous carbon, which makes up most of the dust, is feature- 
less. Both dust species typically condensate at around 2 photo- 
spheric radii at a temperature of around 1500 K (cf. e.g. Ohnaka 
et al., 2007). Molecular shells of C2H2 (bands at 8 - 9 yum and 
1 1 - 14 yum) and HCN (band at 1 1 - 13 jim) at different distances 
and temperatures from the star might be present as well, but are 
not resolved in the low resolution MIDI spectra. The absolute 
flux level of the MIDI spectrum is also reasonable and joins 
with the blackbody curve. In the scheme of Sloan et al. (1998a), 
V Hya is categorized as Red, meaning that a strong continuum 
increase in intensity toward longer wavelengths is superposed on 
the dust features. 

4.3. Visibility modeling results 

The most straight-forward way of interpreting sparsely sam- 
pled interferometric data (visibilities) is by fitting the Fourier 
transform of an assumed brightness distribution of the object. 
Simple size estimations can be obtained from elementary geo- 
metrical models with only a few free parameters. As mentioned 
in the introduction, the definition of a diameter is difficult be- 
cause of its strong wavelength dependence as well as of phase- 
to-phase and long-term cycle-to-cycle variations. On the other 
hand, the size and its dependence on wavelength and pulsation 
phase can give some constraints on the chemical and physical 
mechanisms which are responsible for this appearance, but also 
which layer of the atmosphere is actually observed. This work 



focuses on interpreting the data using simple geometric mod- 
els. In this section, all measurements over all pulsation phases, 
cycles and position angles are taken into account, knowing that 
this potentially biases the interpretation, but will give a first de- 
scription of the size of the structures. In Sect. 5.3, measurements 
are interpreted through each available pulsation phase allowing 
to put constraints on the dynamical mechanisms involved. Self- 
consistent hydrodynamic models will be required for a more re- 
alistic analysis. 

The low surface gravity results in an extended atmosphere 
and temperature structure and the formation of molecular layers 
around late-type stars. Therefore, no sharp transition between 
the star and the circumstellar environment exists. These stars can 
usually be modeled by using a uniform disk (UD) to account for 
the star and a uniform ring to represent a molecular layer (or 
simpler with two uniform disks). However, this would increase 
the numbers of parameters to be fitted, and the later discussed 
pulsation dependence studies (Sect. 5.3) would not be possible. 
Therefore, a simple fully limb-darkened disk (FDD) will be used 
to represent the stellar photosphere and the close molecular envi- 
ronment. One can assume that the fitted FDD diameter will con- 
sequently be located in between the stellar photosphere and the 
outer boundary of the close molecular environment detectable 
with MIDI. 

Uniform disk (UD) model fits are given for comparison 
where appropriate. For some of the stars, a Gaussian function 
is added to constrain the size of the extended dust shell. For the 
stars where this was not possible, i.e. where the dust shell was 
overresolved, only the flux contribution of a surrounding dust 
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shell were determined by not forcing the visibility function to 
be 1 at zero spatial frequency (referred to as relative FDD/UD). 
The total visibility function, V, is written as 

y = I ^FDD ^FDD + (1 - ^FDd) ^Gaussian I , (1) 

with 6fdd the flux contribution of the FDD and the complex 
visibilities (Zhao-Geisler, 2010). For an overresolved dust shell, 

^Gaussian IS Cqual tO ZCrO. 

The final calibrated visibilities are shown in Fig. 3 as func- 
tion of wavelength, whereas the assigned visibility errors are 
given in Tables A. 1 to A.4. As described in paper I, the visibility 
uncertainties are assumed to be the same within each wavelength 
bin to compensate for difficulties in determining them and to en- 
sure a better fitting. The best model fits are displayed in the left 
hand panels of Fig. 4 for each star. The fit and the data, color- 
coded by visual light phase, are plotted as function of spatial fre- 
quency (projected baseline divided by wavelength) for a repre- 
sentative wavelength. It is important to note that diff'erent spatial 
frequencies probe diff'erent regions, i.e. high spatial frequencies 
are sensitive to small regions and low spatial frequencies to ex- 
tended structures. 

The best model parameters with errors are derived by per- 
forming the Levenberg-Marquardt least-squares minimization 
procedure programmed for the interactive data language IDL as 
MPFIT by C. B. Markwardt^^ as described in paper I. The FDD 
diameter as function of wavelength is shown in the right hand 
panels of Fig. 4 for the full data set and selected pulsation phase 
bins (details on the dynamic behavior are discussed in Sect. 5.3). 
The according numerical values are listed in Tables A. 1 to A.4. 
The relative flux contribution of the FDDs and the dust shell 
sizes as function of wavelength are plotted in Fig. 5. The results 
for each star are described in the following in more detail and 
are summarized in Table 5. 

4.3.1. R Aql and R Aqr 

The individual visibility measurements of R Aql and R Aqr 
(Figs. 4a-Fc) have a relatively high scatter. This is caused mainly 
by the low N-band flux of the target and the calibrator and there- 
fore a low signal-to-noise ratio of the fringe signal, but also due 
to the fact that the plot contains observations made at diff'erent 
pulsation phases and pulsation cycles. Intrinsic asynmietries of 
the objects might have an influence on the scatter as well. 

The distribution of the visibilities suggests that two compo- 
nents contribute to the total visibility, since for both stars the val- 
ues are significantly below one at low spatial frequencies. This 
can be attributed to the presence of a dust shell, which was al- 
ready inferred from the MIDI spectra and its comparison with 
the ISO spectra. However, from Figs. 4a-hc, it is clear that the 
size of the dust shell can only be estimated for R Aqr. Due to 
the lack of measurements at very low spatial frequencies this is 
not possible for R Aql (a lower limit is obtained from the limited 
FoV, cf. Sect. 4.2). 

Even though the first zero is not present, the clear decline 
with increasing spatial frequency makes it possible to fit a circu- 
lar FDD or circular UD as second component to the data. A UD 
is in general best suited for a first diameter estimation of a star 
and is therefore often used in the literature. However, the inves- 
tigation of the O-rich stars R Hya and W Hya, presented below, 
suggests that an equally simple geometrical FDD describes the 
observed visibilities in O-rich stars much better in the second 



http://cow.physics.wisc.edu/~craigm/idl/idl.html 



lobe. However, both models are not distinguishable in the case 
of R Aql and R Aqr (Figs. 4a-Fc). 

The behavior of the FDD diameter, ^dd, as function of 
wavelength in the N-band is diff'erent for both stars (Figs. 4b-i-d 
and Tables A.l-hA.2). While for R Aql, ^fdd only sHghtly de- 
creases from (23.1 ± 1.5) mas to (21.3 ± 2.0) mas, going 
from 8 jim to 10 jjm, the relative decrease for R Aqr, from 
(27.4 ± 2.5) mas to (20.9 ± 3.5) mas, is much higher. The 
following relative increase, ^i2//m/^i0//m, beyond 10 jim is how- 
ever similar, being (31 ± 15)% for R Aql and (24 ± 34)% for 
R Aqr. This apparent diameter increase from about 21 mas to 
(27.8 ± 2.0) mas and about 21 mas to (26.0 ± 5.5) mas is equiva- 
lent to an increase from 4.6 AU to 6.2 AU and 5.3 AU to 6.5 AU 
at the assumed distance of R Aql and R Aqr, respectively. 

In comparison, the relative flux contribution of the FDD 
component, 6fdd, decreases from (0.99 ± 0.02) to (0.80 ± 0.02) 
and from (0.57 ± 0.03) to (0.39 ± 0.03) for R Aql and R Aqr, re- 
spectively, going to longer wavelengths (Fig. 5a). This decrease 
reflects the increased flux contribution from the colder surround- 
ing dust shell. In contrast to R Aql, the relative mid-IR flux 
contribution of the star/molecule layer of R Aqr is considerably 
lower because of the larger amount of circumstellar dust bound 
in the symbiotic system. 

For R Aqr, measurements at spatial frequencies lower than 
10 arcsec"^ have larger visibilities than expected for a pure FDD. 
This is illustrated in Fig 4c by fitting a FDD to the measurements 
excluding the ones shortward of 10 arcsec"^ This gives the pos- 
sibility to constrain the inner extension of the silicate dust shell 
interferometrically. However, the Gaussian fit is dominated by 
one point at around 6 arcsec"\ and therefore the result should 
be taken with caution. On the other hand, this single measure- 
ment was taken close in time to the points at around 7.5 arcsec"^ 
and has a position angle not too far off. 

The fitted Gaussian FWHM diameter for R Aqr steadily in- 
creases from (78 ±12) mas to (104 + 7) mas between 8 and 
10 yum, and stays around this value between 10 and 12 //m 
(Fig. 5b). This is equivalent to an increase from 20 AU to 26 AU 
and sets the characteristic silicate dust shell radius at 10 jum at 
a distance of about 7.6 ±1.5 times the photospheric radius. The 
constant or even slightly declining FWHM at wavelengths long- 
ward of 10 Jim could be a hint that the dust shell is truncated 
due to the gravitational conditions set by the orbiting WD with a 
semi-major axis of around 15 AU (60 mas), (cf. Sect. 2.2). With 
the sudden decline of the dust density, cold dust in outer regions, 
probed by wavelengths longer than 10 //m, is less abundant and 
can therefore not be detected. This would be also consistent with 
the low mass-loss rate reported for R Aqr. 

4.3.2. R Hya and W Hya 

The relatively high flux of the target and the calibrator entails 
a comparatively small scatter in the visibility measurements of 
R Hya and W Hya (Fig. 4e, and Fig. 6 in paper 1, respectively). 
However, there is still a considerable spread due to including 
observations obtained at difl'erent pulsations phases, pulsation 
cycles and position angles. From the investigation of W Hya in 
paper I (Fig. 6) and the plot for R Hya in Fig. 4e it can be seen 
that a circular FDD model can well reproduce the visibility am- 
plitude in the second lobe. 

The FDD diameter as function of wavelength in the N-band 
is similar for both stars (Fig. 4f, and Fig. 8 in paper 1). They are 
qualitatively also comparable to R Aql. While the FDD diame- 
ter of R Hya decreases only marginally from (51.6 ± 1.2) mas 
to (46.6 ±0.6) mas between 8 and 10 //m, ^fdd stays almost 
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Fig. 4. Left: Fits of different models to the visibility measurements at a representative wavelength. The visibility data are plotted versus spatial 
frequency and are color-coded by visual light phase. Right: The fitted circular fully limb-darkened disk diameter, ^fdd, as function of wavelength 
for the full data set and selected pulsation phases. The corresponding plots for W Hya can be found in paper I (Fig. 6 and 8, respectively). 
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constant at a value of about (80 ± 1.2) mas for W Hya. At wave- 
lengths longer than 10 yum the apparent diameters gradually in- 
crease again and reach (58.0 ± 0.8) mas and (105 ± 1.2) mas 
at 12 //m, corresponding to a relative increase, ^i2/im/^i0/im, of 
(24 ± 2)% and (31 ± 3)% for R Hya and W Hya, respectively. 
The apparent diameter increase from 47 mas to 58 mas and 
80 mas to 105 mas is equivalent to an increase from 6.1 AU 
to 7.5 AU and 7.1 AU to 9.5 AU at the distance of R Hya and 
W Hya, respectively. 

In order to account for the flux contribution of an extended 
dust shell, the visibility function is not forced to be 1 at zero 
spatial frequency. Figure 5a shows that with longer wavelengths 
the relative flux contribution of the FDD, 6fdd, decreases from 
(0.92 ± 0.03) to about (0.80 ± 0.02) and from (0.85 ± 0.02) to 
about (0.77 ± 0.02) for R Hya and W Hya, respectively, reflect- 
ing the increased flux contribution from the colder surrounding 
dust shell. This is very similar to R Aql. 

4.3.3. V Hya 

The scatter in the visibility measurements in the only C-rich 
AGB star in the sample is acceptable, since the target and the cal- 
ibrator have a comparatively high flux (Fig. 4g). However, there 
is a considerable spread notably at higher spatial frequencies. 
This can probably be attributed to the presence of the temporally 
and spatially changing high velocity outflows (cf. Sect. 2.5). 

A model consisting of two components is necessary to de- 
scribe the visibility measurements. One component is needed 
to explain the moderate decrease of the visibilities from 10 to 
40 arcsec"^ As for the other stars, this could be either a UD or 
an FDD. Both geometrical models cannot be distinguished by 
these observations. As for the 0-rich stars, the FDD will be pre- 
ferred in the following. 

The second component is essential to account for the points 
with high visibilities at spatial frequencies around 7 arcsec"^ 
Similar to 0-rich stars, a Gaussian, representing a carbon- 
rich (AMC -h SiC) dust shell, could be an adequate function. 
However, the sinusoidal visibility variation between 10 and 
20 arcsec"^ suggests a brightness distribution with steep edges. 
A uniform circular ring with a sharp inner and outer edge could 
be a proper function describing a dust shell and causing a si- 
nusoidal modulation. Since the sinusoidal efl'ect is only weak 
at longer wavelengths the dust shell is better represented by a 
Gaussian in that wavelength regime. 



This is shown in Fig. 4g in more detail. At shorter wave- 
lengths the FDD -h ring model fits apparently better (upper panel, 
8.5 yum), while at longer wavelengths the FDD -h Gaussian model 
is better (lower panel, 11.3 yum). However, a comparison of the 
reduced chi square values, xi^ Table A.4 reveals that the 
FDD + Gaussian model is actually the best representation across 
the whole N-band. Nevertheless, the;^^'s are very similar be- 
tween 8 and 9 yum for both models. 

Due to the lack of measurements at very low spatial frequen- 
cies the boundaries of the ring and the FWHM of the Gaussian 
are less well constrained, respectively. Since the fit included 
all measurements, morphological and variability eflTects are ex- 
pected to have a non-negligible influence on the visibility mea- 
surements for this star. In particular, the ring dimensions should 
be taken with caution and the Gaussian FWHM should be used 
as a rough estimation of the inner dust shell boundary. 

The trend of the FDD diameter as function of wavelength 
in the N-band is shown in Fig. 4h. It can be inferred from this 
plot (and Table A.4) that ^fdd decreases from (28.5 ± 2.0) mas 
to (14.9 ± 6.5) mas between 8 and 12 yum with a local maxi- 
mum at around 10.6 yum. This corresponds to a relative decrease 
Ou/um/O^/um of (48 ± 24)%. V Hya is the only star in the sample 
exhibiting a strong diameter decrease from short to long wave- 
lengths. The diameter decrease from 28.5 mas to 14.9 mas is 
equivalent to a decrease from 10.3 AU to 5.4 AU at the distance 
of V Hya. 

The relative flux contribution of the FDD, 6fdd, decreases 
from a maximum of (0.27 ± 0.03) at 8.8 yum to (0.12 ± 0.02) at 
12 yum, reflecting the increased flux contribution from the colder 
surrounding AMC and SiC dust shell detected at longer wave- 
lengths (Fig. 5a and Sect. 5.2). This decrease is similar to the 
0-rich stars in the sample. As for R Aqr, the flux contribution 
of the star/molecule layer is considerably lower. This is proba- 
bly again related to the large amount of dust bound in this sys- 
tem. In addition, the SiC dust shell with a spectral feature in the 
MIDI spectrum at around 1 1.3 yum could also explain part of the 
decrease of the relative flux contribution at longer wavelengths. 

The Gaussian FWHM diameter, Oq, on the other hand 
steadily increases from (99 ± 6) mas to (117 ± 4) mas between 
8 and 12 yum (Fig. 5b). This is equivalent to an increase from 
36 AU to 42 AU and sets the characteristic AMC dust shell ra- 
dius at 10 yum at a distance of 8.7 ± 0.4 times the photospheric 
radius. This dust condensation radius is at a similar distance as 
for the symbiotic 0-rich star R Aqr. 
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The increase of 6q with longer wavelengths is consistent 
with the results for the FDD + ring model. Figure 5b shows that 
the outer ring diameter, ^r,out, increases with increasing wave- 
length as well. Notable is also that the inner ring diameter, in, 
approaches the FDD diameter at longer wavelengths. 

4.4. Wavelength dependence of the diameter 

In the following discussion, one has to keep in mind that the 
compared diameters and diameter ratios (summarized in Table 5) 
are averaged over the position angle, and are averages over the 
pulsation phase or may represent only a certain pulsation phase. 
It is also important to recognize that the FDD diameters obtained 
with MIDI describe a region whose exact location depends on 
the flux contribution of all constituents (continuum photosphere, 
atmospheric molecular layers and nearby dust shells; see next 
section) as function of wavelength and pulsation phase, and not 
only the photosphere of the star. 

The right hand panels of Fig. 6 and Fig. 7 show the obtained 
FDD diameters in relation to interferometric angular diameter 
determinations reported by various authors from the visual to the 
mid-IR (0.6 - 12 jim). They are obtained by fitting a Gaussian, 
a uniform disk or a fully limb-darkened disk to the data. A con- 
version between the models is not performed, since the various 
diameter determinations depend on the number of visibility mea- 
surements and their spatial frequency distribution the authors 
used to fit their model^^. Information on visual phases and posi- 
tion angles (if applicable) of these observations can be found in 
Table B.l (Appendix B). 

The observed apparent diameter changes dramatically within 
the given wavelength range due to the strong wavelength de- 
pendent opacity of the atmospheric constituents (Baschek et al., 
1991; Scholz, 2001). For the 0-rich stars, the measured diame- 
ters in the optical are sensitive to TiO bands. The largest vari- 
ations are around the strongest bands at 712 nm and 670 nm 
with apparent diameter enlargements of up to a factor of two. 
Additionally, light scattered by dust (Norris et al., 2012) might 
also cause a large increase in apparent diameter towards the blue 
(Ireland et al., 2004a). 

In the near-IR, predominantly and CO in diff'erent lay- 
ers are responsible for the wavelength dependence of the diam- 
eter (Hofmann et al., 1998; Jacob et al., 2000). By comparing 
the observations conducted at J, H, K and L band (1.25, 1.65, 
2.16 and 3.8 yum, respectively), one can again infer from the right 



hand panels of Fig. 6 and Fig. 7 that also in the near-IR diameters 
vary by up to a factor of two, and in particular for W Hya that 
there is a complex diameter dependence on pulsation cycle and 
pulsation phase (see Woodruff' et al., 2009, for a more detailed 
discussion on this). This has been reported for a large number of 
AGB stars (e.g. Weiner et al., 2003). Unfortunately, for V Hya, 
only H and K-band observations have been published. Diameter 
determinations for C- stars in these bands are mainly inffuenced 
by the presence of molecular shells consisting of C2H2, HCN, 
CN and CO (Gautschy-Loidl et al, 2004; Paladini et al, 2009). 

Emission from multiple layers and contamination by nearby 
continuum emission have a strong inffuence on determining the 
true photospheric extension of an AGB star through interfero- 
metric measurements. However, a reasonable estimate for the 
photospheric diameter, ^phot, can be obtained from line free mea- 
surements at K-Band. The UD diameter at K-band is approxi- 
mately 1.2 times the true photospheric diameter (cf. e.g. Millan- 
Gabet et al., 2005)^^. The estimated ^phot used in this study 
are listed in Table 5. They were derived by averaging the K- 
band diameters shown in the right hand panels of Fig 6 and 
Fig. 7 (cf. also Table 1 and Table B.l) and dividing them by 1.2 
(^UD,Kband/l-2). The errors are obtained from the standard devia- 
tion. 

Maser observations of diff'erent molecules give additional di- 
ameter constraints. For 0-rich stars, SiO and H2O masers probe 
inner regions where the molecular layers are present and the first 
dust formation takes place, while OH masers trace wind regions 
farther out. Ring diameters for SiO masers, ^sio, were measured 
at 43.1 and 42.8 GHz for R Aql, R Aqr and W Hya by the au- 
thors given in Sect. 2. The averages of these diameters are listed 
in Table 5. A comparison of the location of the SiO masers with 
the photospheric diameters shows that for all three stars the SiO 
masers occur at approximately the same distance from the star. 
The ratios, ^sio/^phot, are between 2.2 and 2.5 (Table 5). 

Except for R Aqr and V Hya, a similar ratio of 2.2 can be 
found if the FDD diameters at 10 yum are compared with the pho- 
tospheric diameters (^FDD,iO/im/^phot)- The SiO masers are there- 
fore almost co-located with the region characterized by the FDD 
diameter (cf. next section). 

For R Aqr and V Hya, the determined FDD diameters lie 
closer to the true photometric diameter. This might be again re- 
lated to the fact that both stars are close binary systems with 
a large dust content and apparently a closer formation of con- 
stituents traced in the N-band. 

In addition, it is notable from Fig. 6d that the apparent di- 
ameters of R Aqr at 3.1 jim and within the L-band are higher 
compared to the K and N-band. This might be caused by an ex- 
tended halo of OH as a result of the ionizing radiation of the 
compact close companion or the existence of hot circumstellar 
dust (Tuthill et al, 2000). 

Millan-Gabet et al. (2003) obtained for the AMC dust 
shell in V Hya in the K-band a FWHM diameter, ^G,2.2/im? 

of 

(35 ± 3) mas (with a ffux contribution of 0.63 ± 0.02). This 
is 3.0 ± 0.3 times smaller compared to the FWHM diameter, 
^G,iO/im, determined at 10 yum, and is therefore consistent with 
the fact that at the N-band colder circumstellar material farther 
out is probed. 

In addition, the AMC dust shell with a FWHM of 
(113 ± 4) mas (and a ffux contribution of 0.85 ± 0.01) deter- 
mined with MIDI for V Hya is in good agreement with the inner 
Gaussian dust shell diameter of (98 ± 2) mas (and a ffux con- 



Empirical conversion factors are approximately: FDD 1.15 UD ^ 
1.68 FWHM. 



However, for C-rich AGB stars, the K-band is probably not a good 
tracer for the true photosphere. 
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tribution of the inner and outer shell of 0.78 ± 0.01) found by 
Townes et al. (2011) at 11.15 fim (narrow bandwidth) with the 
Infrared Spatial Interferometer (IS I). 

It should be also noted that the FDD diameter measured for 
W Hya of (96.3 ± 1.2) mas at around 11.15 jum, tracing the 
AI2O3 dust shell (see next section), is consistent with the av- 
erage UD diameter of (106.2 ± 3.4) mas of the inner dust shell 
determined by Wishnow et al. (2010) with ISI. The inner shell 
flux contribution (including the stellar flux) is with 0.64 ± 0.02 
similar to the value of 0.85 ± 0.01 (Table A.4) measured with 
MIDI as well. 



5. Interpretation and discussion 

Some of the important results from the last sections are summa- 
rized in Table 5 and show that the MIDI observations sample 
the region above the extended pulsating atmosphere where the 
molecular layers are present and probably first seed particles for 
dust formation originate. This study can therefore also help to 
understand the wind acceleration processes in AGB stars. 

The obtained angular diameters as function of wavelength in 
the N-band are significantly diff'erent between the four 0-rich 
stars and the one C-rich star in the sample. R Aql, R Aqr, R Hya 
and W Hya show a moderate to low diameter decrease from 8 
to 10 Jim and a strong diameter increase from 10 to 12 yum, 
while the FDD diameter of V Hya essentially decreases from 
8 to 12 jjm. This is again shown in more detail in Fig. 8a. Since 
diff'erent constituents cause these diff'erences in the shapes, both 
underlying chemistries are discussed separately in the following. 
General conclusions will be derived for oxygen and carbon-rich 
stars concerning their molecular shells and dust production. In 
addition, phase-to-phase variations and asymmetries are investi- 
gated. 

Possible physical mechanism and chemical processes in the 
probed region are constrained by the temperature among other 
things. In particular, the temperature indicates if dust can already 
condensate at the determined FDD distance. A rough estimate 
of the equivalent blackbody temperature can be derived from the 
total flux Sy (in Jy), the flux ratio 6fdd and the diameter ^fdd (in 
rad) measured with MIDI via 

T = — ^ m 

with ly = 4Sy/(l0^^7i0^j^j^) the spectral brightness, v the 
frequency of the observation, h the Planck constant, the 
Boltzmann constant and c the speed of light. The obtained tem- 
peratures are plotted in Fig. 8b for all stars. Due to the use of 
a simple FDD, giving a diameter which lies between the stellar 
photosphere and the outer boundary of the close molecular envi- 
ronment, the real extension of the molecular structure is under- 
evaluated. Therefore, this approximation may overestimate the 
temperature of the close molecular environment. 

R Aql, R Hya and W Hya show very similar trends with tem- 
peratures between 1500 and 2000 K and an error on the order of 
200 K. R Aqr and V Hya exhibit a temperature increase long- 
ward of 1 1 yum, but have considerably higher errors on the order 
of 500 K. The largest error comes from the uncertainty of the to- 
tal flux calibration. However, the temperatures between 10 and 
12 //m, i.e. where a larger FDD diameter were measured, are suf- 
ficiently low for the 0-rich stars so that dust particles are able to 
form (see the following sections). 



5.1. Dust and molecular shells in 0-rich AGB stars 

5.1 .1 . Spectral features and detectability 

The discussion on the spectra in Sect. 4.2 has shown that several 
molecules are present in the upper atmosphere. In the N-band, 
between 8 and 13 //m, strong pure-rotation lines of H2O are 
expected. In addition, SiO exhibits fundamental bands between 
8 and 10 fim (e.g. Decin, 2000). Such quasi-static, warm and 
dense molecular layers close to the star, at typically 2-3 pho- 
tospheric radii (/?phot), have been detected in 0-rich AGB stars 
(e.g. Mennesson et al., 2002; Perrin et al., 2004; Ireland et al., 
2004d; Woodruff et al., 2004; Fedele et al., 2005) and red su- 
pergiant stars (e.g. Perrin et al., 2007). These layers were intro- 
duced earlier to explain spectroscopic observations (e.g. Hinkle 
& Barnes, 1979; Tsuji et al., 1997; Yamamura et al., 1999). 

Besides the classical amorphous silicate feature at 9.7 /im, 
spectral dust features in the N-band probably originate from 
amorphous aluminum oxide (AI2O3) at around 11.5 jim 
(Begemann et al., 1997) and spinel (MgAl204) at around 13 fim 
(e.g. Posch et al., 1999; Fabian et al., 2001). In particular, amor- 
phous AI2O3 provides significant opacity for wavelengths long- 
wards of 10 yum (Koike et al., 1995; Begemann et al., 1997; 
Posch et al., 1999; Woitke, 2006a; Ireland & Scholz, 2006; 
Robinson & Maldoni, 2010). In addition, features of crystalline 
aluminum oxide (corundum, (2-AI2O3) at 12.7 yum and several 
modifications of titanium oxide (e.g. rutile, Ti02) might be 
present as weU (e.g. Posch et al., 1999, 2002). 

It is expected that in 0-rich atmospheres oxides condensate 
first since silicon in silicates have a comparable low electron 
affinity for oxygen (cf. Stencel et al., 1990) and have therefore 
in general a lower condensation temperature due to the lower 
binding energy. The condensation sequence probably starts with 
Ti02 before AI2O3, MgAl204 and Mg^jFcggO forms (Tielens, 
1990; Gail & Sedlmayr, 1998; Jeong et al., 1999; Blommaert et 
al., 2006; Lebzelter et al., 2006; Verhoelst et al., 2009). If most 
of the available high electron affinity metals (primarily Al, Mg, 
Fe, Ca) have been oxidized, silicates, like forsterite (Mg2Si04) 
or olivine (MgQ gFe^ 2^104), may condense and will be absorbed 
onto the seed particles. However, iron-free silicates, nearly trans- 
parent at the stellar ffux maximum at around 1 jim, can survive at 
higher temperatures, thus also condensing at small radii (Norris 
et al., 2012), and are eventually able to trigger the mass-loss 
through photon scattering (Hofner, 2008). Iron-rich silicates, 
which are not transparent to the stellar radiation, thus absorb- 
ing the radiation and heat up, can only exist at larger radii from 
the star. 

Concerning the detectability of the spectral features, two 
characteristics have to be considered. First, if the mass loss is 
relatively low the abundance of silicates will remain low (since 
most of the oxygen is bound in oxides) and the N-band will not 
be dominated by the classical silicate feature (Sogawa & Kozasa, 
1997). This is also related to the geometrical thickness of the 
dust shell as described in Egan & Sloan (2001) (referred to as sil- 
icate dust sequence): spectra showing the classical narrow 10 jim 
feature arise from optically thin shells dominated by amorphous 
silicates while spectra with broad low-contrast emission peaking 
around 11-12 jjm arise from optically and geometrically thin 
shells composed primarily of alumina dust. This might be also 
evidence for different evolutionary stages of 0-rich AGB stars 
(Stencel et al., 1990; Posch et al., 2002; Lebzelter et al., 2006). 
In earlier evolutionary stages the dust mass-loss rates are low 
and aluminum oxide dust dominates while at later stages, when 
eff'ective mass loss has set in, iron-rich silicates will form in large 
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Table 5. Summary of the results. 



Target: 


RAql 


R Aqr 


RHya 


WHya 


VHya 


JjC&L 111UU.C1. 

6fdd layer^ 

^FDD layer^ (mas): 

dust^''' (mas): 


FDD 

0.77 + 0.02 

21.3 + 2.0 
(4.7 AU) 

> 1000 
(> 220 AU) 


FDD + 

Gaussian 

0.35 + 0.03 

20.9 + 3.5 
(5.2 AU) 

105 + 8 
(26.3 AU) 


iV^lClLl V V./ 

FDD 

0.84 + 0.01 

46.6 + 0.6 
(6.1 AU) 

> 1000 
(> 130 AU) 


FDD 

0.78 + 0.01 

79.0+1.2 
(7.1 AU) 

> 1000 
(> 90 AU) 


FDD + 

Gaussian 

0.23 + 0.03 

19.3 + 3.5 
(6.9 AU) 

105 + 4 
(37.8 AU) 


assumed ^^j^^^ (mas): 
assumed 6^.^ (mas): 


y.o + l.Z 

24.2 + 2.7 


1 O Q 1 9 C 
IJ.O + Z.J 

31.7 + 1.2 


ZU.4 + l.U 


Jj. 1 ± L.I 

113 ± 10.5 


1 o 1 ^ f\ 1 
IZ.l + U.J 


^FDD,12//m/ ^FDD,10A/m- 
^Sio/ ^phot- 
t^FDDJO^im/ t'phot- 
lOum / ^„v,^*' 

vj,iw/xiii/ phot 


(31 + 15)% 
2.5 + 0.4 

Z.Z ± U.J 

> 100 


(24 + 34)% 
2.3 + 0.4 

7.6 + 1.5 


(24 + 2)% 

Z.J ± U.l 

>49 


(31 +3)% 
2.2 ± 0.3 
Z.Z ± u.z 
>28 


(18 + 5)%^ 
1 _u n ^ 

l.O + U.J 

8.7 + 0.4 


Phase-to-phase variations: 

^FDD,max/ ^FDD,min- 
^FDD,max / ^FDD,min • 


(18 + 12)% 
(13 + 10)% 


(14 + 19)% 
(18 + 23)% 


(8 + 3)% 
(10 + 10)% 


(13 + 4)% 
(7 + 10)% 




Cycle-to-cycle variations: 






(5 ± 3)% 


(5 ± 4)% 




Asymmetry: 


no 


yes 


no 


yes 


yes 



Notes. ^ At 10 ^m. ^ Lower limits of the dust shell diameter (not Gaussian) for R Aql, R Hya and W Hya are obtained from FoV restrictions 
(cf. Sect. 4.2). ^ Assumed photospheric diameter derived from the K-band diameter (Table 1) divided by 1.2 (cf. e.g. Millan-Gabet et al., 2005). 
^ Assumed SiO ring diameter (average of the values given in Sect. 2, but see also Sect. 4.4). ^ For V Hya, ^G,i2//m/^G,8//m is given. 



amounts farther away from the star and will dominate the dust 
spectra (e.g. Woitke, 2006a). 

The second point concerns the fact that an interferometer al- 
ways acts as a spatial filter. The measured visibility at a given 
baseline, and therefore the determined diameter, is a function 
of the flux contribution (or opacity) of all emitting components 
at that spatial frequency. In addition, only flux contribution of 
components within the field of view (FoV) are measurable. 



5.1 .2. The molecular water shell and close AI2O3 dust shell 

The observations made in this work, in particular the trend of 
the apparent FDD diameter as function of wavelength, summa- 
rized in Fig. 8a, can be understood qualitatively. If they are ad- 
ditionally compared with the results of Ohnaka et al. (2005) and 
Wittkowski et al. (2007) for RR Sco and S Ori, respectively, as 
described in detail for W Hya in paper I, the following conclu- 
sions can be drawn: The overall larger diameter in the mid-IR 
originates from a warm molecular layer of H2O, and the apparent 
gradual increase longward of 10 /im arises most likely from the 
presence of very close AI2O3 dust emitting in this wavelength 
range (labeled in Fig. 8a). 

However, also any other kind of dust compound thermally 
stable in the very close environment of the star and emitting in 
the mid-IR could be responsible for this. The iron-rich silicate 
dust emission does not have any influence since its emitting re- 
gion is over-resolved at the relevant baselines. The comparison 
with the models for RR Sco (Ohnaka et al., 2005) and S Ori 
(Wittkowski et al., 2007) suggests that the partially resolved 
molecular layers are optically thick and that the nearby AI2O3 



dust shell is optically and geometrically thin (cf. also Egan & 
Sloan, 2001). 

The formation of AI2O3 dust at these short distances from 
the stellar surface would be consistent with the empirical results 
by e.g. Lorenz-Martins & Pompeia (2000). The interpretation of 
the apparent diameter increase beyond 10 jim as a result of the 
presence of an AI2O3 dust shell above the molecular water vapor 
shell is supported by the fact that in R Aql and R Hya the spectra 
are characterized by broad oxygen-rich dust emission, meaning 
that oxides are the dominant dust species. In addition, the deter- 
mined equivalent blackbody temperature (Fig. 8b) indicates as 
well that the temperatures are low enough that dust composed of 
titanium and aluminum oxides is stable in the probed region but 
iron-rich silicate dust not. 

For R Aqr, AI2O3 grains are probably present as well, but 
their emission is smeared out by the large amount of amorphous 
silicate dust in the MIDI FoV mainly emitting in the N-band. 
Extinction eff'ects are probably the reason why the diameter in- 
crease longward of 10 jim is first strong but reaches a constant 
value beyond 1 1 //m. In addition, the ionizing radiation of the 
compact companion prevents the formation of large amounts of 
water vapor, meaning that H2O cannot sufficiently form from 
OH and H in the reaction equilibrium. Therefore, the FDD di- 
ameter at around 10 iim is only marginally larger than the pho- 
tospheric diameter. 

5.1 .3. The contribution of SiO 

The only 0-rich star in the sample which does not fit very well 
into the previous interpretation is R Aqr. An overabundance of 
SiO (Angeloni et al., 2007) and a strong silicate dust emis- 
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- TC2H2+ HCN (C-rich) C2H2+ HCN - 




wavelength I (/^m) wavelength I (wm) 

Fig. 8. Left: Normalized FDD diameters of the AGB stars studied in this work (from the right hand panels of Fig. 4). The diameters are normalized 
to one at 10 Errors are approximately the same for all wavelengths and are given on the side. Right: Equivalent blackbody temperature of all 
five AGB stars as function of wavelength (cf. Eq. 2). 



sion could already be inferred from the spectrum and can be 
attributed to the gravitational attraction of the WD in this sym- 
biotic system. The high abundance of SiO molecules and the 
gravitational attraction leads to an eflTective formation of silicate 
dust close to the star and results in a detectable silicate emission 
feature in the MIDI spectrum and a high flux contribution of this 
dust shell. 

At shorter wavelengths, the FDD diameter increases due to 
the presence and increased flux contribution of a SiO molecu- 
lar shell emitting in the 8 to 9 yum range. This is similar to the 
apparent enlargement for the RSG star a Orionis in this wave- 
length regime (Perrin et al., 2007). Even the other 0-rich stars 
show a small diameter enlargement shortwards of 10 jjm due to a 
close SiO molecular shell. A possible relevance of SiO at shorter 
wavelengths is supported by the occurrence of SiO masers in the 
region of the water vapor layer and inner boundary of the pu- 
tative AI2O3 dust shell. Since specific physical conditions are 
necessary in order to exhibit maser emission and these condi- 
tions are solely present at a certain distance from the star, the 
SiO maser shell in R Aqr is also located at a similar distance 
from the star as for the other 0-rich stars in the sample. 

5.1 .4. The outer silicate dust shell 

Since, except for R Aqr, the silicate emission feature is not visi- 
ble in the MIDI spectra, but seen in the ISO spectra, it can be as- 
sumed that there is a surrounding dust shell consistent of mainly 
the more abundant classical silicates, and that even outside the 
FoV of MIDI dust formation is still ongoing. From visibility 
modeling and the FoV some limits on the size of the silicate 
dust shells were obtained. These lower limits are summarized in 
Table 5 and show that the amorphous iron-rich silicate dust is 
located fairly far away from the star. 

However, this could be related for R Aql and R Hya to their 
recent thermal pulse with a short period of enhanced mass loss. 
Such an outer detached shell may contribute significantly to the 
total silicate dust emission. In addition, W Hya is known to have 
a very extended dust shell (Hawkins, 1990). Since the mass-loss 
rates for these three stars are comparably low, the dust shells 
are not very prominent and contribute only weakly to the total 
flux in the mid-IR. In contrast, the emission of the silicate dust 
in the symbiotic system R Aqr contributes substantially to the 



total flux and the characteristic dust condensation radius could 
be determined to approximately 7.6 Rphot- 

5.1 .5. On the wind formation and mass-loss mechanism 

The details of the pulsation-enhanced dust-driven wind mass- 
loss mechanism in 0-rich AGB stars are still under investiga- 
tion (Hofner et al., 2003; Woitke, 2006a; Hofner, 2008; Norris et 
al., 2012). It is believed that the wind is initialized by absorbing 
the momentum of the outward-directed stellar radiation by dust 
and re-emitting it in all directions. The required increase of the 
scale height of the atmosphere is naturally done by the pulsation 
and the formation of shock fronts, accelerating the gas to reach 
the location where it is able to condense into dust grains. Since 
amorphous AI2O3 has only a moderate abundance, it can proba- 
bly not be responsible to initiate the mass-loss (cf. e.g. Woitke, 
2006a). AI2O3 can exist close to the star without inducing mass 
loss. Unfortunately, the observations presented here cannot give 
a conclusive answer which dust species are responsible for trig- 
gering a high mass loss. 

As described in paper I, the scattering oflf large iron-free sil- 
icate grains, as proposed by Hofner (2008), cannot be verified 
either, since they were not detected. However, this might be re- 
lated to their low abundance making their emission smeared out 
by H2O, SiO, and AI2O3. The low abundance of iron-free sili- 
cates is not critical for driving the wind since a small fraction of 
Si condensed into Mg2Si04 would be already suflftcient (Hofner, 
2008). Recent observations in scattered light in the near-IR by 
Norris et al. (2012) support the existence of large iron-free sil- 
icates at around 2 stellar radii. Since for this first study no ra- 
diative transfer modeling associated to self-consistent dynamic 
modeling was done, no strong conclusions about the grain mix- 
ture and wind formation in the close environment can be drawn. 

Not many other wind acceleration mechanisms come in 
mind. There might be the possibility that small amounts of car- 
bon grains play a role as speculated by Hofner & Andersen 
(2007), but with the side eflfect that this would lead to a too 
high infrared excess. Despite its very low abundance, Ti02 might 
have some relevance as well (Posch et al., 2002), but is not de- 
tectable with these observations. In this context, it might be in- 
teresting to know if scattering on AI2O3 grains is important. Also 
the role of large amounts of water vapor in molecular shells and 
the radiation pressure on water molecules may need more de- 
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tailed calculations. In addition, metallic Fe might have some ef- 
fects as well (McDonald et al., 2010). 

It is clear that quantitative modeling is necessary to support 
the above findings, in particular, if the derived constituents of the 
close molecular and dust layers (H2O, SiO, AI2O3) can really 
provide sufficient opacities to explain the observed diameter de- 
pendence on wavelength, in particular if AI2O3 could cause the 
apparent diameter increase beyond 10 fim. In order to quantify 
the results, it will be necessary to apply dynamic atmospheric 
models of e.g. Ireland & Scholz (Ireland & Scholz, 2006; Ireland 
et al., 2011) or Hofner et al. (see Lebzelter et al., 2010; Sacuto 
et al., 2011). Even if AI2O3 dust is not consistently included in 
these models so far, and has often to be added ad hoc, this will 
give more detailed insight into the physical processes at work 
there. 

5.2. Dust and molecular shells in C-rich AGB stars 

5.2.1. Spectral features 

The study of the oxygen-rich stars in the sample has shown 
that dense molecular layers, extending to up to 2 photospheric 
radii, are present in these stars in agreement with findings in 
other 0-rich AGB stars. In contrast, the physical properties of 
the outer atmosphere of carbon stars and their temporal varia- 
tions have not yet been probed well. Hydrostatic and dynamic 
model atmospheres of carbon-rich Miras as well as non-Mira 
carbon stars fail to explain observed spectra longward of about 
5 fim (J0rgensen et al., 2000; Gautschy-Loidl et al., 2004). These 
models predict very strong absorption due to C2H2 and HCN at 
7 and 14 fim, but the observed ISO spectra of carbon stars show 
only weak absorption due to these molecular species. However, 
Aoki et al. (1998) and Aoki et al. (1999) proposed that emis- 
sion from extended warm molecular layers, containing C2H2 and 
HCN, may be responsible for the discrepancy between the ISO 
spectra and the models. Both species are among the most abun- 
dant in C-rich atmospheres. 

While in the K and H band the spectrum is contaminated 
with CO bands, absorption features of C2H2 and HCN, as well 
as emission from SiC at 11.3 yum (Kozasa et al., 1996) and fea- 
tureless amorphous carbon (AMC) dust, play an important role 
in the N-band. In addition, molecules such as silane (11 jim) and 
ammonia (10.7 //m) may also provide significant opacities in C- 
rich stars (Monnier et al., 2000, and references therein). 

The only observations with MIDI of classical C-rich stars 
were conducted by Ohnaka et al. (2007) for V Oph and Sacuto 
et al. (2011) for R Scl. Ohnaka et al. (2007) modeled V Oph 
with a dust shell, consisting of AMC (85%) and SiC (15%), with 
a multi-dimensional Monte Carlo code (Ohnaka et al., 2006), 
and added a polyatomic molecular layer, consisting of C2H2 
and HCN, where the opacities of these molecules were calcu- 
lated with an appropriate band model assuming LTE. Sacuto et 
al. (2011) applied a self-consistent dynamic model atmosphere 
to their data (Hofner et al., 2003; Mattsson et al., 2008, 2010). 
With this model, molecular shells of C2H2 and HCN above the 
stellar photosphere could be probed. However, the dusty envi- 
ronment could not be adequately sampled and a complementary 
model was applied (a hydrostatic stellar model plus a dust radia- 
tive transfer code). 

5.2.2. The dust and molecular shells 

The FWHM of the Gaussian distribution in V Hya (Fig. 5b) 
slightly increases by about 20% longwards of 10 fim probably 



due to the presence of SiC dust surrounding the star and emit- 
ting in the N-band. The FDD diameter as function of wavelength 
(Fig. 4h) shows a strong increase of 30% shortward of 9 fim 
which could be related to the presence of molecular layers of 
C2H2 and HCN (see Ohnaka et al., 2007; Sacuto et al., 2011). 
These molecular layers would then lie at about 2 photospheric 
radii. Noteworthy, the equivalent blackbody temperature short- 
ward of 9 fim gets low enough 2000 K, Fig. 8b), supporting 
the fact that molecules are able to form. 

There is, however, a discrepancy of the FDD diameter at 
wavelengths longward of 1 1 //m. With the presence of the men- 
tioned molecular shells one would also expect some increase of 
the FDD diameter in this wavelength part. There are possibly 
several reasons for this discrepancy. First of all, the uncertain- 
ties in fitting a FDD (determined by the visibilities with spa- 
tial frequencies larger than 10 arcsec"^) are much higher in this 
wavelength regime due to a much higher scatter of the visibilities 
(cf. Fig. 4h and Table A.4). 

This scatter comes from combining observations done at dif- 
ferent epochs and position angles. For example, asymmetries are 
very pronounced in V Hya (cf. Sect. 5.4). The position angles 
where the observations were made are similar to the position an- 
gles of the high velocity wind (cf. Fig. 6g and e.g. Hirano et al., 
2004) and temporal fluctuations of the density and emissivity of 
the massive outflows can be expected. 

Therefore, it can only be concluded that SiC and certainly 
AMC dust is present at the very close circumstellar environ- 
ment. The dust-driven wind formation scenario in C-rich stars is 
to a certain amount satisfactorily solved (e.g. Gail & Sedlmayr, 
1987; Hofner et al, 2003; Gautschy-Loidl et al., 2004; Nowotny 
et al., 2005; Woitke, 2006b; Wachter et al., 2008). Radiation 
pressure on carbon grains is efficient enough to initiate the wind. 

The presence of C2H2 and HCN in V Hya cannot be con- 
firmed with certainty due to measurement uncertainties and ap- 
plying only a simple geometric model. However, also the visi- 
bility as function of wavelength (Fig 3d) is very similar to the 
ones of V Oph and R Scl (Fig. 1 in Ohnaka et al. (2007) and 
Fig. 3 Sacuto et al. (2011), respectively), suggesting the same 
atmospheric conditions and the presence of those layers. 

V Hya shows some similarities to R Aqr except for the un- 
derlying chemistry. As for R Aqr, the average N-band FDD di- 
ameter is only about 1 .6 times larger than the photospheric diam- 
eter, and the dust condensation radius is at about 8 times the pho- 
tospheric radius. However, it is not expected that the K-band di- 
ameter is a good tracer for the true photospheric radius in C-rich 
stars, and the condensation of AMC dust could occur at smaller 
radii, i.e. at higher temperatures. 

As for the oxygen-rich stars, it is necessary to use dynamic 
model atmospheres (e.g. Sacuto et al., 2011), computed with 
stellar parameters corresponding to V Hya, for a better under- 
standing of the physical processes responsible for the molecule 
and dust formation close to the star. 



5.3. Dynamic behavior of the close environment 

As a result of the intrinsic pulsation, mid-IR flux variations on 
the order of 20%, with a phase offset of about 0.15, were de- 
tected in this study except for R Aql (Sect. 4.1). This section 
investigates how this relates to changes of the apparent angular 
sizes and flux contributions of the close stellar atmosphere, and 
how this influences or triggers the dust formation and mass loss. 
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visual light phase visual light phase 



Fig. 9. Phase-to-phase variations of all oxygen-rich stars in the sample. The left penal shows the normalized FDD diameter as function of visual 
light phase and the right panel the normalized relative flux contribution of the FDD as function of visual light phase (normalized in respect to the 
full data set). It was not possible to obtain reliable results for the C-rich star V Hya. 



Table 6. Phase-to-Phase Variations. 



Phase Phase No of Norm.^ FDD Norm.^ Flux 
Bin^ Range Obs. Diameter Contribution 



RAql: 

1 0.875-0.125 8 

2 0.125-0.375 9 0.95 + 0.03 0.98 + 0.06 

3 0.375-0.625 9 1.08 + 0.08 1.10 + 0.06 

4 0.625-0.875 6 1.11 + 0.09 1.01 + 0.06 



R Aqr: 

1 0.750-0.250 12 0.93 + 0.10 0.92 + 0.12 

2 0.250-0.750 14 1.06 + 0.10 1.08 + 0.12 



RHya: 



1 


0.000- 


-0.200 









2 


0.200- 


-0.400 


14 


1.04 + 0.02 


0.94 + 0.06 


3 


0.400- 


-0.600 


27 


1.01 + 0.02 


1.04 + 0.06 


4 


0.600- 


-0.800 


20 


0.96 + 0.01 


1.00 + 0.06 


5 


0.800- 


-0.000 


3 


1.00 + 0.04 


1.00 + 0.06 



WHya: 

1 0.875-0.125 23 1.09 + 0.02 1.03 + 0.06 

2 0.125-0.375 42 0.96 + 0.02 0.96 + 0.06 

3 0.375-0.625 10 1.02 + 0.04 1.06 + 0.06 

4 0.625-0.875 - 



Notes. ^ Cf. right hand panels of Fig. 1 . ^ Normalized in respect to the 
full data set. ^ No reliable fit could be obtained. 



5.3.1 . Size and flux variations 

Cycle-to-cycle variations were studied for all oxygen-rich stars 
for the phase ranges given in the left hand panels of Fig. 1. For 
this, the best geometrical model was fitted to the visibilities with 
the relative fluxes fixed to the value obtained for the full data set 
(cf. Fig. 5a). No cycle-to-cycle investigation could be done for 
V Hya due to the fact that observations are not repeated at the 
same phase in consecutive cycles. 

However, reliable cycle-to-cycle variations could be ob- 
tained for R Hya and W Hya (for W Hya see paper I). Both 
stars show a maximum variation of the N-band averaged and 
normalized FDD diameter, ^fdd (normalized in respect to the 
full data set), on the order of (5 ± 4)%, thus much lower than 
phase-to-phase variations as shown below. This verifies that the 



folding of consecutive pulsation cycles into one cycle is an ac- 
ceptable assumption for the following phase-to-phase analysis. 
For R Aql and R Aqr, the uncertainties are higher and no reli- 
able trend could be found. 

For studying the phase-to-phase behavior of the diameter and 
the relative flux contribution of the FDD, the pulsation phase 
of each star is divided into several phase bins as shown in the 
right hand panels of Fig. 1 and listed in Table 6. As for the 
cycle-to-cycle investigation, the best geometrical model is then 
fitted to the data within each bin. For the fit to R Aqr, the 
Gaussian FWHM is fixed to the value obtained for the full data 
set (cf. Fig. 5b). No reliable fits could be obtained for V Hya be- 
cause of a very unfavorable spatial frequency coverage. In sev- 
eral phase bins the FDD diameter could not been fitted or the 
uncertainties were too high due to potential asymmetries and un- 
known cycle-to-cycle variations (cf. Sect. 5.4) even after reduc- 
ing the number of phase bins. Therefore, V Hya is not considered 
in the following discussion. 

For a simple comparison between phase bins, the obtained 
FDD diameters, ^fdd (cf. right hand penals of Fig. 4), and rela- 
tive flux contributions, 6fdd, are averaged over the N-band wave- 
length regime (8-12 jim). Since the chromatic trend of ^fdd and 
6fdd does not change considerably between phases, and in gen- 
eral other error sources are more prominent (see below), this ap- 
proach can be appropriate for a first investigation. The similar 
chromatic trend found for each subsequent phase bin indicates 
also that the probed diflTerent layers behave similarly as function 
of pulsation phase. 

The N-band averaged diameters, 0, and relative flux contri- 
butions, 6, are plotted as function of visual phase in Fig. 9, and 
are listed in Table 6 in comparison to the full data set values. The 
errors are in general relatively high because of the low number 
of visibility points used for these fits. In addition, the large phase 
binning reduces the true amplitude of the variation somewhat. It 
should be also noted that for R Aqr only two phase bins are used 
and that there are no fit results and observations available for 
R Aql and R Hya at visual maximum, respectively. 

Except for R Hya, the FDD diameter is largest at or af- 
ter the visual minimum (Fig. 9a and Table 5). The increase, 
^FDD,max/^FDD,min, bctwccn minimum and maximum is on the or- 
der of 15%. While the layer representing the photosphere in a 
Mira variable is in the mid-infrared largest near the stellar lumi- 
nosity maximum (Weiner et al., 2003), the more outer pulsating 
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layers reach their maximum extension with an increasing phase 
delay. This is in agreement with the movement of a mass-shell at 
around 2 stellar radii as modeled by Nowotny et al. (2010) (their 
Fig. 1 and 2) for a carbon star. However, for R Hya a different 
trend was found making this interpretation uncertain. The large 
error bars also reflecting this uncertainty. 

In contrast, the variation of the relative flux contribution as 
function of visual light phase shows a more consistent behavior. 
The flux contribution is lower at around visual maximum than 
at around visual minimum (Fig. 9b and Table 5). The increase, 
^FDD,max/^FDD,min5 bctwccn minimum and maximum is on the or- 
der of 10 to 15%, but also with large uncertainties and therefore 
not strongly significant. 

The location of the outer silicate dust shell as function of 
visual phase was studied for the symbiotic system R Aqr. For this 
the Gaussian FWHM was not fixed to the full data set value. At 
visual minimum, the averaged Gaussian FWHM, Oq, was larger 
than at visual maximum by about (24 ± 19)%. 

5.3.2. Dust and wind formation 

Taking the above findings at face value, one infers that new dust 
forms around or after the visual minimum when the pulsating 
layer at around 2 photospheric radii reaches its maximum exten- 
sion. Since a higher relative flux contribution of the FDD around 
the visual minimum is approximately canceled out by the lower 
absolute flux, the temperature in this region will be lower too. 

Calculating the equivalent blackbody temperature according 
to Eq. 2, the temperatures around the visual minimum are on 
average (150 ± 100) K lower and at around the visual maximum 
on average (150+ 100) K higher than the temperatures averaged 
over the cycles shown in the right panel of Fig. 8. This means that 
in particular in the 11-12 fim range, where the putative AI2O3 
dust shell is traced, the temperatures around the visual minimum 
are getting as low as the condensation temperature of AI2O3 of 
approximately 1400 K, showing that new AI2O3 dust could form. 

The phase behavior would also be consistent with the theo- 
retical predictions of Ireland & Scholz (2006). In their model, 
the largest rate of grain growth occurs between phases 0.3 and 
0.7, similar to what is predicted for carbon-rich stars (cf. Fig. 2 
in Nowotny et al., 2010). 

While the observations made in this work show that AI2O3 
dust might newly form around the visual minimum, it cannot 
be determined conclusively how this relates to the acceleration 
of the wind and the mass-loss rate. As mentioned before and in 
paper I, AI2O3 can exist in the upper atmosphere without being 
a wind-driver. No conclusion can be given concerning the rel- 
evance of micrometer- sized iron- free silicates for the wind for- 
mation since this needs the use of radiative transfer modeling as- 
sociated to self-consistent dynamic modeling taking opacities of 
different realistic chemical compounds into account. However, it 
can be speculated that the relevant constituents that are respon- 
sible for the wind acceleration are formed along with AI2O3. 

5.4. Asymmetries 

The good uv-coverage for R Hya and W Hya, shown in the left 
hand panels of Fig. 6 and Fig. 1 in paper I, makes it possible to 
directly fit a elliptical FDD to the data. For R Aql, R Aqr and 
V Hya, the uv-planes are not filled well. Departure from spher- 
ical symmetry was here studied by simply comparing visibility 
measurements at different position angles taken closely in time 



and with similar projected baselines. Differential phases were 
not used but will be investigated in more detail in future studies. 

As described in detail in paper I, W Hya shows a small elon- 
gation in the mid-IR. The position angle (FA) and axis ratio of an 
elliptical FDD were determined to (11 ± 20)° and 0.87 ± 0.07, 
respectively. The asynmietry was explained by a possible en- 
hanced dust concentration along an N-S axis. As for W Hya, 
an elliptical FDD was fitted to the full data set of R Hya, but 
also to subsets of the full data set to test for time dependencies. 
However, no departure from spherical symmetry could be re- 
vealed within the measurement uncertainties. Any intrinsic de- 
viation from sphericity is less then 10%. 

While for W Hya contradictory position angles were found 
in the literature, the non-detection of an asymmetry in R Hya is 
consistent with the results of Ireland et al. (2004a) and Monnier 
et al. (2004). Both authors did not find any departures from sym- 
metry within the measurement errors in the optical and near-IR 
(680-940 nm and 2.26 fim, respectively). 

Because both stars have similar intrinsic properties but a 
different asynmietric appearance, they are excellent targets for 
studying the origin and formation of asymmetries in more de- 
tail, e.g. if they are related to different evolutionary phases. By 
comparing both stars, it could be investigated why W Hya shows 
an asymmetry and R Hya not, and in particular why in W Hya 
asymmetries appear different at different scales and how they are 
connected. 

Figure 10 shows selected visibility measurements for R Aql 
and V Hya to investigate their asymmetry. In order to do so, 
groups of points with similar spatial frequencies are compared 
with the best model fit to the full data set serving as reference^^. 
In R Aql, it can be seen that the visibilities measured at vi- 
sual light phase 0.17 with a PA of 75° (crosses) lie on the full 
data set fit as well as basically also the measurements taken at 
an almost perpendicular PA of 129° (diamonds). Therefore, it 
can be concluded that there is no departure from symmetry in 
the close vicinity of the star at around 2.2 photospheric radii 
within the measurement uncertainties. In contrast, it can be in- 
ferred that at a distinct different visual light phase the visibilities 
(squares) are smaller and therefore probing a region appearing 
larger (cf. Sect. 5.3.1). 

In V Hya an opposite effect can be seen. Points with a PA 
around 65° (squares) have a higher visibility than expected for 
a circular FDD and points with PAs from 130 to 180° a lower 
visibility. This seems not be related to the pulsation since at dif- 
ferent light phases but similar PAs the visibilities are similar. 
This means that a strong asymmetry already sets in in the very 
close stellar environment at less then 2 photospheric radii with 
an elongation along a North-South direction, therefore, a simi- 
lar orientation as found by e.g. Tsuji et al. (1988); Kahane et al. 
(1988, 1996); Sahai et al. (2003) and particular also by Lagadec 
et al. (2005) in the mid-infrared (cf. Sec. 2.5). In contrast to 
R Aql, where the symmetry is seen over the whole N-band 
wavelength regime, the asymmetry in V Hya appears stronger 
at longer wavelengths (10 to 12 fim). 

For R Aqr, only one set of visibilities could be compared. In 
the shorter wavelength regime (8 to 10 yc/m), the visibilities taken 
at 2007-10-07 05:54 at a PA of 79° (phase 0.87, Bp^oj = 40.3 m) 



Note that a comparably large phase and spatial frequency range is 
used. This allows a more reliable analysis due to comparing groups of 
points instead of only single measurements with large errors. However, 
it can not be excluded that this is biased by the fact that different spa- 
tial frequencies probe different parts of the envelope. Therefore, only a 
qualitative analysis is given here. 
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Fig. 10. Selected visibility measurements to investigate the asymmetries in R Aql (left) and V Hya (right). The visual light phase and position 
angle are given in parenthesis for each data point. See text for more details. 



are larger than the ones taken at 2007-10-07 01:16 (phase 0.87, 
Bproj = 41.3 m) at a PA of 59° (cf. Table C.2). This is a hint that 
the distribution of SiO strongly deviates from spherical symme- 
try. From this investigation it follows that three out of the five 
stars observed show asymmetries. 

6. Summary 

Five AGB stars, namely R Aql, R Aqr, R Hya, W Hya, and 
V Hya, were observed over a period of nearly three years cov- 
ering several pulsation cycles. All stars were for the first time 
monitored in the mid-IR (8-13 jim) with the interferometric in- 
strument MIDI. These high-resolution observations are sensitive 
to the structure of the stellar atmosphere, consisting of the con- 
tinuum photosphere and overlying molecular layers, as well as 
to the properties of the dust shell. 

The results and conclusions can be summarized as follows: 

- Except for R Aql, a periodic flux variation is found in the 
N-band. The typical semi-amplitudes are on the order of 
20% to 30% for the Miras and 10% for the SR/Mira variable 
W Hya. The amplitude is probably also related to the dust 
content in the system. The mid-IR maximum occurs always 
after the visual maximum at visual phase 0.15 ± 0.05. 

- In the MIDI spectrum of R Aqr, the silicate dust emission 
feature could be clearly identified, while it was not present in 
the spectra of the other 0-rich stars. This has been attributed 
to the field-of-view of MIDI smaller than the one of ISO. The 
expected SiC dust emission feature is detected in the MIDI 
spectrum of the carbon star V Hya. 

- The obtained visibility data were best fitted by a fully limb- 
darkened disk (FDD), representing the central star and the 
close atmospheric layers, and a spherical Gaussian distribu- 
tion, representing the extended dusty environment. 

- The trend of the FDD diameter as function of wavelength is 
similar for all 0-rich stars in the sample. The apparent size 
is almost constant between 8 and 10 yum (21.3 ± 2.0 mas, 
20.9 ± 3.5 mas, 46.6 ± 0.6 mas, and 79.0 ± 1.2 mas, for 
R Aql, R Aqr, R Hya, and W Hya, respectively) and grad- 
ually increases at wavelengths longer than 10 yum. The en- 
largement is for all stars between 25% and 30%. In contrast, 
the relative flux contribution of the FDD decreases, reflect- 
ing the increased flux contribution from the surrounding dust 
sheU. 

- The measured FDD diameters at 10 yum are about 2.2 times 
larger than the photospheric diameters (except for R Aqr and 



V Hya). This ratio is consistent with observations of other 
AGB stars. The smaller ratio of about 1.6 for R Aqr and 

V Hya can be understood from the gravitational interaction 
with a compact companion. 

- The obtained diameter trends for the 0-rich stars are 
compared with the findings of Ohnaka et al. (2005) and 
Wittkowski et al. (2007) for RR Sco and S Ori, respectively. 
It can be concluded that the overall larger diameter in the 
mid-IR originates from a warm molecular layer of H2O, and 
the gradual increase long ward of 10 yum can be most likely 
attributed to the contribution of a spatially resolved close 
AI2O3 dust shell. A contribution of SiO shortward of 10 jim 
can be seen as well, maybe except for W Hya. 

- These observations of a larger sample of stars than available 
before confirm previous results, and emphasize the need for 
dynamic models able to treat the formation of dust in a self- 
consistent way. In particular, the close dust shell composed 
of AI2O3 well below the distance at which the iron-rich sili- 
cate dust shell is traced, could now be revealed for a signifi- 
cant number of stars. 

- A comparison with available SiO maser ring diameters 
shows that the SiO maser emission is co-located with the 
region probed with MIDI. 

- For the carbon rich star V Hya it can be concluded that SiC 
and AMC dust is present in the very close circumstellar en- 
vironment. The chromatic shape of the Gaussian FWHM 
shows an increase around the SiC feature (11.3 yum). The 
apparent FDD diameter as function of wavelength shows a 
strong increase of more than 30% shortward of 9 yum, which 
could be a hint that close molecular layers consisting of 
C2H2 and HCN at about 2 photospheric radii are probed. 
However, a simultaneous increase at wavelengths longward 
of 1 1 yum could not been detected, but this might be related 
to large fitting uncertainties. 

- It is found that the iron-rich silicate dust shells are located 
fairly far away from the star, at radii larger than 30 times 
the photospheric radius. Only the silicate dust shell of R Aqr 
is located closer to the star at around 8 photospheric radii, 
which is most probably related to the interaction with a close 
companion. 

- Phase-to-phase variations of the FDD diameter and the rel- 
ative flux contribution of the 0-rich stars are on the order 
of 15% and 10-15%, respectively. It was found that the rela- 
tive flux contribution was highest at around visual minimum, 
while there was no clear behavior detectable for the FDD di- 
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ameter as function of visual light phase. Due to the large 
uncertainties this is not strongly significant and no firm con- 
clusions can be drawn concerning the mass-loss mechanism. 
Nevertheless, one can speculate that new dust forms around 
visual minimum. 

- Maximal diameter changes due to cycle-to-cycle variations 
are on the order of 6%, and are therefore lower than phase- 
to-phase variations. This could only be determined for R Hya 
and W Hya. 

- The close environment geometry could be investigated for 
all stars. R Aql and R Hya do not show any deviation from 
sphericity, while R Aqr, W Hya and V Hya do. 

Most of the proposed goals could be accomplished. 
However, a few items remain open. First, no clear connection 
could be established between the dust condensation and wind 
formation. Second, the size of the iron-rich silicate dust shell 
could not be satisfactorily constrained for 3 of the 4 oxygen-rich 
stars due to the lack of observations at short baselines. Third, the 
pulsation phase coverage allowed only for some of the stars an 
analysis of the apparent diameter with time. 

The majority of the above conclusions rely on fits of a rela- 
tively simple FDD to the observational data. However, this can 
only be an approximation, even if this model represents the data 
quite well. In order to obtain a more meaningful diameter and the 
exact composition, it is therefore necessary to compare with hy- 
drodynamic model atmospheres of e.g. Ireland & Scholz (Ireland 
& Scholz, 2006; Ireland et al., 2011) and Hofner et al. (see 
Lebzelter et al., 2010; Sacuto et al, 2011). 

Another big reward can be gained if these results can be 
combined with kinematic information. Complementary high- 
resolution radio and mm observations, tracing the relevant con- 
stituents and sampHng the dust formation zone probed with 
MIDI, are therefore necessary. 
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Appendix A: Visibility model fit results for R Aql, R Aqr, R Hya and V Hya 

Table A.l. R Aql: circular fully limb-darkened disk (FDD) and circular uniform disk (UD). 
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Table A.2. R Aqr: circular fully limb-darkened disk (FDD) + circular Gaussian (G) and circular uniform disk (UD) + circular Gaussian (G). 
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Notes. ^ mean visibility error used for the corresponding wavelength bin (cf. Sect. 4.3); reduced chi square 
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Table A.3. R Hya: circular fully limb-darkened disk (FDD) and circular uniform disk (UD). 



A 
(jim) 



AV^ 



circular FDD 

^FDD Flux 6fdD 

(mas) 



circular UD 

^UD Flux 6uD 

(mas) 



8.12 


0.012 


51.6 


+ 


1.2 


0.92 


+ 


0.03 


13 


45.6 


+ 


1.2 


0.91 


+ 


0.02 


19 


8.33 


0.012 


50.5 


+ 


0.7 


0.92 


+ 


0.02 


13 


43.4 


+ 


1.2 


0.89 


+ 


0.03 


18 


8.54 


0.012 


50.4 


+ 


0.8 


0.92 


+ 


0.02 


14 


43.2 


+ 


0.8 


0.89 


+ 


0.02 


19 


8.71 


0.011 


49.9 


+ 


0.7 


0.94 


+ 


0.01 


13 


42.9 


+ 


0.6 


0.91 


+ 


0.02 


18 


8.87 


0.011 


48.6 


+ 


0.8 


0.93 


+ 


0.02 


12 


42.1 


+ 


0.6 


0.91 


+ 


0.02 


17 


9.07 


0.012 


47.7 


+ 


0.5 


0.92 


+ 


0.01 


10 


41.1 


+ 


0.6 


0.90 


+ 


0.01 


13 


9.26 


0.011 


47.0 


+ 


0.5 


0.90 


+ 


0.01 


9 


40.5 


+ 


0.6 


0.88 


+ 


0.01 


11 


9.45 


0.013 


46.4 


+ 


0.6 


0.88 


+ 


0.01 


7 


40.3 


+ 


0.6 


0.86 


+ 


0.02 


9 


9.63 


0.011 


46.8 


+ 


0.6 


0.86 


+ 


0.01 


9 


40.3 


+ 


0.6 


0.84 


+ 


0.01 


11 


9.78 


0.011 


46.7 


+ 


0.5 


0.85 


+ 


0.01 


9 


40.3 


+ 


0.6 


0.83 


+ 


0.01 


11 


9.92 


0.011 


46.6 


+ 


0.6 


0.84 


+ 


0.01 


8 


40.3 


+ 


0.6 


0.82 


+ 


0.01 


10 


10.09 


0.011 


46.1 


+ 


0.5 


0.82 


+ 


0.01 


8 


40.3 


+ 


0.6 


0.81 


+ 


0.01 


10 


10.26 


0.011 


46.8 


+ 


0.5 


0.82 


+ 


0.01 


9 


40.9 


+ 


0.8 


0.81 


+ 


0.01 


11 


10.42 


0.009 


48.8 


+ 


0.8 


0.81 


+ 


0.01 


14 


42.2 


+ 


0.6 


0.80 


+ 


0.01 


18 


10.58 


0.009 


50.0 


+ 


0.6 


0.81 


+ 


0.01 


16 


43.2 


± 


0.6 


0.79 


+ 


0.01 


21 


10.71 


0.010 


50.5 


+ 


0.7 


0.80 


+ 


0.01 


16 


43.7 


+ 


0.6 


0.79 


+ 


0.01 


20 


10.84 


0.009 


51.3 


+ 


0.8 


0.80 


+ 


0.01 


19 


44.3 


+ 


0.6 


0.78 


+ 


0.01 


24 


10.99 


0.008 


52.8 


+ 


0.7 


0.79 


+ 


0.01 


21 


45.5 


+ 


0.6 


0.77 


+ 


0.01 


28 


11.14 


0.008 


53.4 


+ 


0.7 


0.78 


+ 


0.01 


21 


46.0 


+ 


0.8 


0.77 


+ 


0.02 


27 


11.29 


0.010 


54.1 


+ 


0.7 


0.79 


+ 


0.01 


16 


46.6 


+ 


0.8 


0.77 


+ 


0.01 


21 


11.43 


0.011 


54.4 


+ 


0.8 


0.79 


+ 


0.01 


13 


46.9 


+ 


0.8 


0.77 


+ 


0.01 


16 


11.55 


0.011 


55.5 


+ 


0.8 


0.80 


+ 


0.01 


14 


47.8 


+ 


0.8 


0.78 


+ 


0.01 


17 


11.67 


0.011 


56.3 


+ 


0.8 


0.80 


+ 


0.02 


15 


48.5 


+ 


1.0 


0.78 


+ 


0.02 


19 


11.81 


0.008 


57.3 


+ 


0.8 


0.80 


+ 


0.01 


29 


49.3 


+ 


1.2 


0.78 


+ 


0.02 


37 


11.95 


0.007 


58.0 


+ 


0.8 


0.81 


+ 


0.02 


40 


49.9 


+ 


1.0 


0.79 


+ 


0.02 


51 



Table A.4. V Hya: circular fully limb-darkened disk (FDD) + circular Gaussian (G) and circular fully limb-darkened disk (FDD) + circular 
uniform ring (R). 



A 
(//m) 



AV^ 



circular FDD -i- circular Gaussian 
Flux 6fdd 



CPDD 

(mas) 



(mas) 



Xi 



2,b 



CPDD 

(mas) 



circular FDD -i- circular ring 

Flux 6fdd ^R,in ^R,out 

(mas) (mas) 



2,b 



8.12 


0.011 


28.5 


+ 


2.0 


0.22 


+ 


0.03 


99 + 6 


8 


35.3 


+ 


2.0 


0.32 


+ 


0.04 


51 


+ 


4 


124 


+ 


6 


14 


8.33 


0.011 


26.5 


+ 


2.0 


0.25 


+ 


0.03 


100 


+ 7 


8 


31.9 


+ 


2.0 


0.33 


+ 


0.04 


49 


+ 


4 


125 


+ 


6 


12 


8.54 


0.012 


24.5 


+ 


2.0 


0.26 


+ 


0.03 


101 


+ 5 


8 


30.0 


+ 


2.0 


0.34 


+ 


0.04 


50 


+ 


4 


127 


+ 


6 


12 


8.71 


0.011 


22.9 


+ 


2.0 


0.27 


+ 


0.01 


101 


+ 5 


9 


28.7 


+ 


2.0 


0.35 


+ 


0.04 


51 


+ 


4 


129 


+ 


6 


12 


8.87 


0.011 


21.8 


+ 


2.0 


0.27 


+ 


0.01 


102 


+ 5 


10 


27.9 


+ 


2.0 


0.35 


+ 


0.04 


51 


+ 


4 


131 


+ 


6 


12 


9.07 


0.011 


21.0 


+ 


2.0 


0.27 


+ 


0.03 


102 


+ 5 


11 


27.7 


+ 


2.0 


0.35 


+ 


0.04 


51 


+ 


4 


133 


+ 


6 


15 


9.26 


0.010 


20.4 


+ 


2.5 


0.26 


+ 


0.03 


102 


+ 4 


13 


27.7 


+ 


2.0 


0.34 


+ 


0.04 


51 


+ 


4 


135 


+ 


6 


17 


9.45 


0.010 


20.2 


+ 


2.5 


0.25 


+ 


0.01 


103 


+ 5 


14 


27.9 


+ 


2.0 


0.34 


+ 


0.04 


51 


+ 


4 


138 


+ 


6 


20 


9.63 


0.009 


19.8 


+ 


3.0 


0.24 


+ 


0.03 


105 


+ 5 


18 


28.1 


+ 


2.0 


0.33 


+ 


0.04 


52 


+ 


4 


140 


+ 


6 


26 


9.78 


0.009 


19.2 


+ 


3.0 


0.23 


+ 


0.03 


105 


+ 5 


19 


28.1 


+ 


2.0 


0.32 


+ 


0.04 


51 


+ 


4 


143 


+ 


6 


29 


9.92 


0.008 


19.3 


+ 


3.5 


0.23 


+ 


0.03 


105 


+ 4 


20 


28.6 


+ 


2.0 


0.32 


+ 


0.04 


51 


+ 


4 


145 


+ 


6 


30 


10.09 


0.008 


19.8 


+ 


3.5 


0.22 


+ 


0.01 


105 


+ 4 


25 


29.3 


+ 


2.0 


0.31 


+ 


0.04 


50 


+ 


4 


147 


+ 


6 


39 


10.26 


0.007 


20.4 


+ 


3.5 


0.21 


+ 


0.01 


106 


+ 4 


26 


30.4 


+ 


2.0 


0.30 


+ 


0.04 


50 


+ 


4 


149 


+ 


6 


41 


10.42 


0.006 


20.6 


+ 


3.5 


0.20 


+ 


0.03 


107 


+ 4 


33 


31.0 


+ 


2.0 


0.29 


+ 


0.04 


49 


+ 


4 


152 


+ 


6 


56 


10.58 


0.006 


20.6 


+ 


3.5 


0.18 


± 


0.03 


108 


+ 5 


35 


31.0 


+ 


2.0 


0.27 


+ 


0.04 


46 


+ 


4 


155 


+ 


6 


65 


10.71 


0.007 


20.6 


+ 


4.0 


0.17 


+ 


0.03 


109 


±5 


29 


31.1 


+ 


2.0 


0.26 


+ 


0.04 


44 


+ 


4 


157 


+ 


6 


57 


10.84 


0.006 


19.8 


+ 


4.5 


0.17 


+ 


0.03 


111 


+ 4 


36 


30.7 


+ 


2.0 


0.25 


+ 


0.04 


43 


+ 


4 


160 


+ 


6 


73 


10.99 


0.005 


19.2 


+ 


4.5 


0.16 


+ 


0.01 


112 


+ 4 


41 


30.4 


+ 


2.0 


0.24 


+ 


0.04 


41 


+ 


4 


162 


+ 


6 


85 


11.14 


0.006 


18.1 


+ 


5.5 


0.15 


+ 


0.01 


113 


+ 4 


37 


29.2 


+ 


2.0 


0.22 


+ 


0.04 


37 


+ 


4 


165 


+ 


6 


77 


11.29 


0.006 


17.5 


+ 


6.0 


0.15 


+ 


0.03 


114 


+ 3 


33 


28.6 


+ 


2.0 


0.21 


+ 


0.04 


34 


+ 


4 


167 


+ 


6 


70 


11.43 


0.006 


16.9 


+ 


5.5 


0.14 


+ 


0.01 


115 


+ 4 


33 


27.7 


+ 


2.0 


0.20 


+ 


0.04 


29 


+ 


4 


169 


+ 


6 


71 


11.55 


0.006 


16.3 


+ 


5.0 


0.14 


+ 


0.03 


115 


+ 5 


33 


27.8 


+ 


2.0 


0.19 


+ 


0.04 


28 


+ 


4 


171 


+ 


6 


73 


11.67 


0.006 


16.8 


+ 


5.5 


0.13 


+ 


0.01 


117 


+ 4 


30 


28.8 


+ 


2.0 


0.19 


+ 


0.04 


28 


+ 


4 


173 


+ 


6 


69 


11.81 


0.006 


15.0 


+ 


6.5 


0.12 


+ 


0.03 


117 


+ 4 


34 


28.9 


+ 


2.0 


0.18 


+ 


0.04 


28 


+ 


4 


175 


+ 


6 


79 


11.95 


0.007 


14.9 


+ 


6.5 


0.12 


+ 


0.01 


117 


+ 4 


30 


25.5 


+ 


2.0 


0.13 


+ 


0.04 


64 


+ 


4 


146 


+ 


6 


70 



Notes. ^ mean visibility error used for the corresponding wavelength bin (cf. Sect. 4.3); reduced chi square 
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Appendix B: Diameter details 



Table B.l. Details of the diameter measurements presented in the right hand panels of Fig. 6 and Fig. 7. 



R Aql: Author 


Model 


Visual phase 


Comments 


A: Haniffetal. (1995) 


UD-kG 


0.06 


UD and Gaussian at 700 and 710 nm 


B: Millan-Gabet et al. (2005) 


UD 


0.90 


UD diameter at H and K band 


C: Ragland et al. (2006) 


UD 


0.70 


UD diameter at H band 


D: Hofmann et al. (2000) 


UD 


0.17 


UD diameter at K band 


E: van Belle et al. (1996) 


UD 


0.90 & 0.31 


upper for phase 0.31 & lower for phase 0.90 (K band) 


F: this work 


FDD 


0.0-1.0 


full data set with error bars 


R Aqr: Author 


Model 


Visual phase 


Comments 


A: Tuthill et al. (1999) 


UD+G 




UD and Gaussian at 830 nm 


B:Tuthill etal. (2000) 


UD 


0.12 & 0.68 


phase 0.68 for upper/lower point at 2.2/3.1 /im 


C: Millan-Gabet et al. (2005) 


UD 


0.40 


UD diameter at J, H and K band 


D: Ragland et al. (2006) 


UD 


0.30 


UD diameter at H band 


E: Ragland et al. (2008) 


UD 


0.60-1.11 


UD diameters around the H band at different phases 


F: Mennesson et al. (2002) 


UD 


0.41 & 0.51 


UD diameter at K and L band 


G: van Belle et al. (1996) 


UD 


0.34 & 0.57 


upper for phase 0.34 & lower for phase 0.57 (K band) 


H: this work 


FDD 


0.0-1.0 


full data set with error bars 




iVlULlCl 


visual piid&c 


l^UllllllCllLh 


A: Ireland et al. (2004a) 


G 


0.62 


upper for PA 123° & lower for FA 72° 


B: Haniffetal. (1995) 


UD-kG 


0.28 


UD and Gaussian at 700, 833 and 902 nm 


C: Monnier et al. (2004) 


UD 


0.50 


UD diameter at K band 


D: Millan-Gabet et al. (2005) 


UD 


0.80 


UD diameters at J, H and K band 


E: this work 


FDD 


0.0-1.0 


full data set with error bars 


W Hya: Author 


Model 


Visual phase 


Comments 


A: Lattanzi et al. (1997) 


UD 


0.64 


two DPrnpnrliriilar axp<s (n x 86 PA 14^° SS'^ nm^ 


B: Ireland etal. (2004a) 


G 


0.44 


upper curve for PA 120° and lower curve for PA 252° 


C: Haniffetal. (1995) 


UD 


0.04 


UD diameter at 700 and 710 nm 


D: Tuthill et al. (1999) 


G 


0.04 


elliptical Gaussian at 700 and 710 nm (77 ^ 0.94, PA ^ 93°) 


E: Monnier et al. (2004) 


UD 


0.50 


UD reported as a bad fit to the data (K band) 


F: Woodruff etal. (2009) 


UD 


0.58-1.53 


curves for phase 0.58 (middle), 0.79 (lower) and 1.53 (upper) 


G: Woodruff etal. (2008) 


UD 


0.50-1.00 


mean of multiple measurements in given phase range (J, H, K, L) 


H: Millan-Gabet et al. (2005) 


UD 


0.60 


UD diameter at H and K band 


1: Wishnow et al. (2010) 


UD 


0.1 - 1.1 


inner UD dust shell diameter at 11.15 fim 


J: Bedding etal. (1997) 


UD 


0.43 


UD diameter at 1.45 fim 


K: this work 


FDD 


0.0-1.0 


full data set with error bars 


V Hya: Author 


Model 


Visual phase 


Comments 


A: Ragland et al. (2006) 


UD 


0.10 


UD diameter at H band 


B: van Belle etal. (1999) 


UD 


0.00 


UD diameter at H band 


C: Millan-Gabet et al. (2003) 


UD 




UD diameter at K band 


D: this work 


FDD 


0.0-1.0 


full data set with error bars 



Appendix C: Observation log 

Table C.l. AT stations used in this study (cf. http://www.eso.org/sci/facilities/paranal/telescopes/vlti/configuration/index.html). 



Name 


AT stations 


Ground length 


Ground PA 


Resolution (at 10 fim) 


A 


EO-GO 


16.006 m 


71.020° 


128 mas 


B 


GO-HO 


31.998 m 


70.998° 


64 mas 


B* 


AO-DO 


32.011m 


71.014° 


64 mas 


C 


EO-HO 


48.004 m 


71.005° 


42 mas 


D 


DO-HO 


64.005 m 


71.012° 


32 mas 


E 


DO-Gl 


71.557 m 


134.443° 


28 mas 


F 


HO-Gl 


71.555 m 


7.576° 


28 mas 
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Table C.2. Detailed observation logs of R Aql, R Aqr, R Hya, and V Hya. 



Observation log of R Aql: 



Date 


UT time 


pa 


AT'' 


Disp'^ 


Cal^ 


B^ (m) 


PA^ (°) 


Phases 




QF 


2007-04-17 


09:06:23.779 


P79 


B 


GRISM 


2 of 3 


29.05 


74.77 


0.78 


53.87 


used 


2007-04-25 


08:09:12.338 


P79 


D 


GRISM 


2 of 2 


54.90 


74.83 


0.81 


50.84 


used 


2007-05-09 


09:51:15.000 


P79 


D 


GRISM 


of 1 


63.39 


71.22 


0.85 


52.97 


not used 


2007-06-19 


08:25:08.000 


P79 


B 


GRISM 


of 2 


28.98 


66.54 


0.02 


41.70 


not used 


2007-06-20 


03:28:40.000 


P79 


c 


GRISM 


4 of 6 


34.14 


73.90 


0.02 


41.23 


used 


2007-06-20 

A^\J\J 1 \J\J A^\J 


04:42:25.000 


P79 


B 


GRISM 


of 1 


28.64 


74.81 


0.02 


52.51 


not used 


2007-06-21 


08:30:51.000 


P79 


A 


GRISM 


of 2 


14.04 


65.19 


0.02 


39.22 


not used 


2007-06-22 


08-05-24 000 


P79 


c 


GRISM 


1 of 1 


44.27 


67.37 


0.02 


43.09 


used 


2007-10-02 


23:37:29.000 


P80 


A 


PRISM 


1 of 1 


16.01 


72.65 


0.40 


56.10 


used 


2007-10-05 


00:24:43.000 


P79 


B 


GRISM 


1 of 2 


31.24 


70.22 


0.41 


50.57 


used 


2007-10-05 


23:33:31.000 


P80 


B 


PRISM 


1 of 1 


32.00 


72.56 


0.41 


55.85 


used 


2007-10-06 


23:43:58.000 


P80 


c 


PRISM 


2 of 2 


47.87 


71.90 


0.42 


54.73 


used 


2007-10-07 


00:11:23.000 


P80 


c 


PRISM 


2 of 2 


47.08 


70.51 


0.42 


51.29 


used 


2008-03-06 


09:16:57.330 


P80 


D 


PRISM 


1 of 1 


31.21 


69.28 


0.98 


28.48 


used 


2008-03-06 


09:26:46.124 


P80 


D 


PRISM 


1 of 1 


33.42 


70.36 


0.98 


30.49 


used 


2008-03-25 


09-31-47 779 


P79 


A 


GRISM 


of 1 


12.33 


74.45 


0.05 


45.38 


not used 


2008-04-02 


09:00:39.945 


P81 


A 


PRISM 


3 of 3 


12.34 


74.46 


0.07 


45.44 


used 


2008-04-02 


09-39-51 946 


P81 


c 


PRISM 


3 of 4 


41.50 


74.83 


0.07 


51.25 


used 


2008-04-04 


09:30:40.946 


P81 


B 


PRISM 


1 of 1 


27.57 


74.82 


0.08 


51.08 


used 


2008-04-04 


09-39-46 960 


P81 


B 


PRISM 


1 of 1 


28.17 


74.82 


0.08 


52.22 


used 


2008-04-27 


08:24:58.000 


P81 


D 


PRISM 


1 of 1 


58.66 


74.76 


0.17 


53.96 


used 


2008-04-27 


08:36:46.960 


P81 


D 


PRISM 


1 of 1 


59.45 


74.69 


0.17 


55.04 


used 


2008-05-23 


05:54:09.000 


P79 


D 


GRISM 


1 of 1 


52.03 


74.69 


0.26 


47.67 


used 


2008-05-26 


05:17:58.000 


P81 


D 


PRISM 


1 of 1 


49.44 


74.46 


0.27 


43.68 


used 


2008-06-01 


06:08:05.138 


P81 


A 


PRISM 


4 of 4 


70.28 


128.8 


0.30 


54.02 


used 


2008-06-01 

A^\J\J\J \J\J \J A. 


06-16-17 960 

% • A. 1 %y\j\j 


P81 


A 


PRISM 


4 of 4 


69.88 


128.8 


0.30 


54.79 


used 


2008-06-01 


06:56:41.946 


P81 


A 


PRISM 


4 of 4 


67.14 


129.2 


0.30 


57.02 


used 


2008-06-01 


07-04-33 959 


P81 


A 


PRISM 


of 4 


66.46 


129.4 


0.30 


57.13 


not used 


2008-06-06 


07:46:19.000 


P81 


A 


PRISM 


of 2 


15.96 


71.90 


0.32 


54.21 


not used 


2008-06-06 


09:40:04.000 


P81 


c 


PRISM 


Oof 1 


41.28 


64.31 


0.32 


36.72 


not used 


2008-06-07 


06:16:17.945 


P81 


C 


PRISM 


1 of 1 


45.96 


74.42 


0.32 


56.43 


used 


2008-06-07 


06:39:05.000 


P81 


B 


PRISM 


3 of 3 


31.67 


73.66 


0.32 


57.11 


used 


2009-04-20 


08:01:43.946 


P83 


B 


PRISM 


3 of 3 


25.57 


74.61 


0.49 


47.18 


used 


2009-04-23 


06:40:13.946 


P83 


C 


PRISM 


2 of 6 


28.53 


72.10 


0.50 


34.71 


used 


2009-04-23 


07:53:25.946 


P83 


C 


PRISM 


6 of 6 


38.78 


74.66 


0.50 


47.72 


used 


2009-05-02 


09:38:45.000 


P83 


D 


PRISM 


3 of 3 


63.97 


72.88 


0.53 


56.23 


used 


2009-06-03 


05:35:12.000 


P83 


A 


PRISM 


1 of 3 


13.88 


74.85 


0.65 


51.00 


used 


2009-06-03 


06:52:23.000 


P83 


A 


PRISM 


lof3 


15.69 


74.01 


0.65 


57.07 


used 


2009-06-26 


04:31:29.000 


P83 


D 


PRISM 


1 of 2 


58.70 


74.75 


0.73 


54.09 


used 


2009-07-03 


05:21:54.000 


P83 


B 


PRISM 


1 of 2 


31.88 


73.27 


0.76 


56.82 


used 


Observation log of R Aqr: 

Date UT time 


pa 


AT'' 


Disp'^ 


Cal^ 


B^ (m) 


PAf n 


Phase^ 




QF 



2007- 


■05-16 


10:06:28.000 


P79 


D 


GRISM 


1 of 1 


51.09 


53.81 


0.50 


50.38 


used 


2007- 


■06-18 


09:46:44.000 


P79 


A 


GRISM 


2 of 4 


15.53 


67.42 


0.59 


74.20 


used 


2007- 


■06-19 


09:34:00.000 


P79 


B 


GRISM 


Oof 2 


30.79 


66.67 


0.59 


72.48 


not used 


2007- 


■06-20 


08:25:38.000 


P79 


C 


GRISM 


4 of 6 


42.81 


61.15 


0.59 


58.66 


used 


2007- 


■06-21 


06:59:42.000 


P79 


A 


GRISM 


2 of 2 


11.32 


45.06 


0.59 


40.09 


used 


2007- 


■07-02 


07:18:59.000 


P79 


B 


GRISM 


2 of 3 


26.64 


56.61 


0.62 


54.29 


used 


2007- 


■07-27 


10:07:03.000 


P79 


D 


GRISM 


1 of 1 


57.40 


78.05 


0.68 


60.39 


used 


2007- 


■10-04 


06:26:11.000 


P79 


A 


GRISM 


Oof 1 


12.59 


79.91 


0.86 


49.04 


not used 


2007- 


■10-04 


02:12:15.000 


P80 


B 


PRISM 


2 of 2 


30.08 


64.81 


0.86 


68.14 


used 


2007- 


■10-04 


04:10:49.000 


P80 


B 


PRISM 


2 of 2 


31.80 


73.85 


0.86 


77.57 


used 


2007- 


■10-04 


05:28:51.000 


P80 


A 


PRISM 


1 of 1 


14.50 


77.84 


0.86 


61.91 


used 


2007- 


■10-05 


05:16:14.000 


P79 


B 


GRISM 


Oof 2 


29.57 


77.35 


0.86 


63.83 


not used 


2007- 


■10-06 


01:57:28.000 


P80 


C 


PRISM 


2 of 2 


44.52 


63.84 


0.87 


66.67 


used 


2007- 


■10-07 


05:54:46.000 


P79 


C 


GRISM 


1 of 1 


40.31 


79.13 


0.87 


53.48 


used 


2007- 


■10-07 


01:16:29.000 


P80 


C 


PRISM 


2 of 2 


41.29 


58.76 


0.87 


58.52 


used 


2007- 


■10-07 


02:00:53.165 


P80 


A 


PRISM 


1 of 1 


14.88 


64.02 


0.87 


68.23 


used 


2007- 


■11-30 


00:31:14.000 


P80 


D 


PRISM 


Oof 1 


47.51 


74.32 


0.01 


76.85 


not used 


2007- 


■12-02 


00:23:50.000 


P80 


D 


PRISM 


Oof 3 


63.42 


74.20 


0.01 


76.78 


not used 


2008- 


■05-22 


09:16:47.000 


P81 


D 


PRISM 


Oof 3 


48.10 


49.62 


0.43 


45.12 


not used 


2008- 


■05-30 


07:15:31.945 


P80 


D 


PRISM 


Oof 5 


36.48 


21.26 


0.47 


24.69 


not used 


2008- 


■05-30 


09:49:08.137 


P80 


D 


PRISM 


2 of 5 


55.36 


59.13 


0.47 


59.52 


used 
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Date 


UT time 


pa 




Disp^ 


Car' 


(m) 


PA*' (°) 


Phase^ 






2008-06-06 


08:15:59.945 


P81 


A 


PRISM 


Oof 2 


11.72 


47.73 


0.49 


44.70 


not used 


2008-06-06 


09:53:36.946 


P81 


C 


PRISM 


1 of 1 


44.07 


63.13 


0.49 


66.54 


used 


2008-06-07 


09:21:11.000 


P81 


B 


PRISM 


3 of 3 


28.21 


60.38 


0.49 


60.30 


used 


2008-06-07 


09:54:00.945 


P81 


B 


PRISM 


3 of 3 


29.58 


63.60 


0.49 


67.47 


used 


2008-09-25 


05:26:02.000 


P81 


D 


PRISM 


1 of 1 


61.26 


76.25 


0.78 


69.53 


used 


2008-09-28 


05:26:37.330 


P81 


A 


PRISM 


1 of 1 


15.23 


76.49 


0.78 


66.90 


used 


2008-12-28 


01:14:16.000 


P82 


A 


PRISM 


Oof 1 


12.25 


80.18 


0.01 


43.58 


not used 


2008-12-31 


00:43:23.000 


P82 


B 


PRISM 


1 of 2 


26.09 


79.50 


0.02 


47.72 


used 


2008-12-31 


01:14:25.000 


P82 


B 


PRISM 


1 of 2 


23.39 


80.52 


0.02 


40.95 


used 


2009-01-01 


00:12:38.000 


P82 


C 


PRISM 


1 of 1 


39.63 


79.35 


0.03 


53.72 


used 


2009-05-25 


09:28:26.000 


P83 


D 


PRISM 


Oof 1 


50.56 


53.10 


0.39 


50.22 


not used 


2009-06-03 


09:04:03.000 


P83 


A 


PRISM 


2 of 3 


13.01 


55.02 


0.42 


52.70 


used 


ZUUy-UD-Uj 


nQ-/in-/i^ nnn 

U>'.4U.4j.UUU 


roj 




PPTQA/r 
r KloiVl 


i 01 i 




DU. JO 


U.4Z 




used 


2009-07-03 


06:48:28.000 


P83 


B 


PRISM 


1 of 2 


24.83 


51.87 


0.50 


48.71 


used 


2009-07-30 


08:40:36.000 


P83 


D 


PRISM 


1 of 1 


63.29 


74.41 


0.57 


75.51 


used 


2009-08-14 


08:16:39.000 


P83 


C 


PRISM 


Oof 1 


45.81 


76.36 


0.62 


68.58 


not used 


2009-08-14 


08:38:16.000 


P83 


A 


PRISM 


Oof 1 


14.85 


77.26 


0.62 


63.93 


not used 


2009-08-14 


09:30:22.000 


P83 


B 


PRISM 


Oof 1 


26.11 


79.49 


0.62 


52.30 


not used 


Observation log of R Hya: 




















Date 


U 1 time 


pa 


AT^ 


Disp'^ 


Cal^ 


B^ (m^ 




Phase^ 






2007-04-12 


02:08:09.000 


P79 


A 


GRISM 


2 of 3 


13.48 


48.88 


0.60 


52.64 


used 


2007-04-17 


07:28:33.338 


P79 


B 


GRISM 


2 of 3 


24.73 


86.35 


0.61 


49.77 


used 


2007-04-18 


07:10:38.780 


P79 


A 


GRISM 


3 of 3 


12.92 


85.18 


0.61 


52.93 


used 


2007-04-22 


01:35:49.071 


P79 


D 


GRISM 


8 of 8 


54.75 


50.43 


0.62 


55.11 


used 


2007-04-22 


07:13:32.070 


P79 


D 


GRISM 


8 of 8 


48.03 


87.13 


0.62 


47.83 


used 


2007-04-25 


01:23:03.779 


P79 


E 


GRISM 


Oof 2 


63.50 


123.9 


0.63 


54.00 


not used 


2007-04-25 


04:23:39.780 


P79 


F 


GRISM 


2 of 2 


71.55 


9.98 


0.63 


84.59 


used 


2007-05-16 


23:52:59.000 


P79 


B 


GRISM 


Oof 1 


27.96 


52.53 


0.69 


52.30 


not used 


2007-06-20 


01:20:22.000 


P79 


C 


GRISM 


6 of 6 


46.25 


77.60 


0.78 


76.25 


used 


2007-06-21 


22:54:49.000 


P79 


C 


GRISM 


1 of 1 


46.13 


62.66 


0.78 


71.24 


used 


2007-07-02 


23:17:47.000 


P79 


C 


GRISM 


5 of 5 


47.95 


69.99 


0.81 


86.20 


used 


2007-07-02 


01:42:07.000 


P79 


C 


GRISM 


Oof 5 


41.92 


82.60 


0.81 


60.55 


not used 


2007-07-03 


22:54:48.000 


P79 


A 


GRISM 


Oof 2 


15.91 


68.30 


0.82 


81.96 


not used 


2007-07-03 


00:28:49.000 


P79 


A 


GRISM 


2 of 2 


15.57 


76.81 


0.82 


76.34 


used 


2007-07-04 


23:58:37.000 


P79 


B 


GRISM 


Oof 4 


31.72 


74.45 


0.82 


82.30 


not used 


2007-07-04 


02:19:46.000 


P79 


B 


GRISM 


Oof 4 


24.68 


86.41 


0.82 


50.24 


not used 


2007-12-26 


08:39:14.000 


P79 


D 


GRISM 


1 of 1 


51.86 


44.87 


0.29 


45.73 


used 


2008-02-20 


07:53:42.000 


P80 


D 


PRISM 


4 of 4 


63.89 


69.61 


0.44 


85.31 


used 


2008-02-20 


08:42:11.000 


P80 


D 


PRISM 


6 of 4 


63.50 


74.33 


0.44 


83.27 


used 


2008-02-21 


09:02:26.000 


P79 


D 


GRISM 


Oof 2 


62.50 


76.42 


0.44 


77.80 


used 


2008-02-22 


07:02:15.000 


P80 


B* 


PRISM 


Oof 6 


31.22 


64.71 


0.45 


75.47 


not used 


2008-02-22 


07:54:30.156 


P80 


B* 


PRISM 


Oof 6 


32.00 


70.49 


0.45 


87.16 


not used 


2008-03-02 


06:36:33.000 


P80 


E 


PRISM 


2 of 2 


70.92 


129.9 


0.47 


77.67 


used 


2008-03-02 


07:52:58.000 


P80 


F 


PRISM 


1 of 2 


71.55 


10.27 


0.47 


84.61 


used 


2008-03-02 


08:28:33.000 


P80 


D 


PRISM 


2 of 2 


62.29 


76.75 


0.47 


76.56 


used 


2008-03-13 


04:43:47.000 


P80 


B 


PRISM 


7 of 8 


29.03 


56.32 


0.50 


61.84 


used 


2008-03-13 


07:00:13.945 


P80 


B 


PRISM 


7 of 8 


31.97 


72.36 


0.50 


86.68 


used 


2008-03-13 


09:09:41.945 


P80 


B 


PRISM 


7 of 8 


27.22 


83.52 


0.50 


57.36 


used 


2008-03-14 


05:39:07.000 


P80 


A 


PRISM 


4 of 6 


15.57 


64.36 


0.50 


75.34 


used 


2008-03-14 


07:02:11.000 


P80 


A 


PRISM 


4 of 6 


15.96 


73.21 


0.50 


85.40 


used 


2008-03-14 


08:51:47.000 


P80 


A 


PRISM 


4 of 6 


13.98 


82.60 


0.50 


60.53 


used 


2008-03-25 


06:08:58.000 


P80 


C 


PRISM 


6 of 6 


47.97 


72.15 


0.53 


87.51 


used 


2008-03-25 


07:59:34.945 


P80 


C 


PRISM 


6 of 6 


42.97 


81.63 


0.53 


62.57 


used 


2008-03-31 


04:48:31.000 


P80 


C 


PRISM 


1 of 1 


47.29 


66.29 


0.54 


79.04 


used 


2008-04-01 


03:50:22.946 


P81 


B 


PRISM 


2 of 2 


29.75 


58.88 


0.55 


66.69 


used 


2008-04-02 


06:42:37.945 


P81 


C 


PRISM 


4 of 4 


46.11 


77.84 


0.55 


72.93 


used 


2008-04-02 


08:01:14.946 


P81 


A 


PRISM 


3 of 3 


13.26 


84.39 


0.55 


55.05 


used 


2008-04-03 


07:04:02.946 


P81 


B 


PRISM 


2 of 2 


29.69 


79.96 


0.55 


67.16 


used 


2008-04-28 


01:55:35.946 


P81 


D 


PRISM 


6 of 6 


58.79 


57.60 


0.62 


64.72 


used 


2008-04-28 


02:57:39.946 


P81 


E 


PRISM 


2 of 2 


70.99 


130.1 


0.62 


78.85 


used 


2008-04-28 


04:06:01.945 


P81 


F 


PRISM 


1 of 2 


71.55 


9.69 


0.62 


85.24 


used 


2008-04-28 


07:13:08.945 


P81 


D 


PRISM 


6 of 6 


44.07 


89.17 


0.62 


42.84 


used 


2008-05-25 


00:16:55.000 


P81 


D 


PRISM 


2 of 2 


59.74 


59.31 


0.69 


66.43 


used 


2008-07-03 


01:31:24.000 


P81 


A 


PRISM 


1 of 2 


14.10 


82.30 


0.79 


61.43 


used 


2008-07-03 


02:43:24.000 


P81 


C 


PRISM 


2 of 2 


33.75 


88.69 


0.79 


45.14 


used 


2008-07-06 


03:20:45.000 


P81 


D 


PRISM 


2 of 2 


35.68 


93.62 


0.80 


34.14 


used 


2009-01-16 


08:03:46.946 


P82 


F 


PRISM 


4 of 5 


71.54 


-7.84 


0.32 


57.04 


used 
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Date 


UT time 


pa 




Disp^ 


Car' 


(m) 


PA*' (°) 


Phase^ 






2009-01-16 


08:34:11.000 


P82 


D 


PRISM 


1 of 2 


58.75 


57.54 


0.32 


63.95 


used 


2009-01-17 


07:50:45.946 


P82 


E 


PRISM 


2 of 2 


63.99 


124.0 


0.32 


54.99 


used 


2009-01-21 


07:17:41.000 


P81 


B 


PRISM 


5 of 5 


26.59 


47.49 


0.33 


51.06 


used 


2009-01-21 


09:02:59.000 


P82 


B 


PRISM 


5 of 5 


31.10 


64.15 


0.33 


74.98 


used 


2009-01-22 


07:23:42.945 


P81 


C 


PRISM 


4 of 5 


40.45 


48.91 


0.34 


53.31 


used 


2009-01-22 


07:36:00.946 


P81 


C 


PRISM 


4 of 5 


41.40 


51.25 


0.34 


56.10 


used 


2009-01-25 


08:13:46.945 


P82 


A 


PRISM 


8 of 8 


14.94 


59.33 


0.34 


67.36 


used 


2009-01-25 


08:50:11.945 


P82 


A 


PRISM 


8 of 8 


15.56 


64.18 


0.34 


75.65 


used 


2009-01-26 


09:11:32.000 


P82 


B 


PRISM 


1 of 1 


31.71 


67.44 


0.35 


81.40 


used 


2009-01-27 


08:13:44.945 


P82 


C 


PRISM 


4 of 4 


45.26 


60.42 


0.35 


69.14 


used 


2009-01-27 


08:39:07.946 


P82 


C 


PRISM 


4 of 4 


46.52 


63.77 


0.35 


74.92 


used 


2009-01-27 


08:46:50.123 


P82 


C 


PRISM 


4 of 4 


46.83 


64.72 


0.35 


76.68 


used 


2009-02-16 


07:39:19.000 


P82 


D 


PRISM 


2 of 2 


63.27 


66.95 


0.40 


79.20 


used 


2009-03-16 


08:08:24.000 


P82 


D 


PRISM 


2 of 3 


59.80 


79.60 


0.48 


68.84 


used 


2009-03-16 


08:36:32.000 


P81 


D 


PRISM 


2 of 3 


56.78 


82.00 


0.48 


62.43 


used 


2009-04-20 


07:34:42.945 


P83 


B 


PRISM 


2 of 3 


23.02 


88.14 


0.57 


45.29 


used 


2009-04-23 


05:28:10.000 


P83 


C 


PRISM 


6 of 6 


45.43 


78.85 


0.58 


71.30 


used 


2009-04-23 


07:05:14.946 


P83 


C 


PRISM 


6 of 6 


36.81 


86.57 


0.58 


49.27 


used 


2009-04-24 


06:05:35.000 


P83 


B 


PRISM 


2 of 2 


28.16 


82.31 


0.58 


61.90 


used 


2009-05-02 


01:44:52.000 


P83 


D 


PRISM 


3 of 3 


59.37 


58.63 


0.60 


65.64 


used 


2009-05-03 


01:50:41.946 


P83 


E 


PRISM 


2 of 2 


68.81 


126.7 


0.61 


67.85 


used 


2009-05-03 


06:19:20.000 


P83 


F 


PRISM 


1 of 3 


70.55 


22.78 


0.61 


50.75 


used 


2009-05-24 


23:13:50.000 


P83 


D 


PRISM 


1 of 1 


53.58 


48.26 


0.66 


51.89 


used 


2009-06-04 


23:30:42.000 


P83 


C 


PRISM 


3 of 3 


44.27 


58.02 


0.69 


64.64 


used 


2009-06-04 


03:46:00.000 


P83 


A 


PRISM 


1 of 2 


13.39 


84.09 


0.69 


56.99 


used 



Observation log of V Hya: 



Date 


UT time 


pa 


AT^ 


Disp'^ 


Cal^ 


B^ (m) 


PA^ (°) 


Phase^ 




QP 


2007-04-12 


03:45:42.000 


P79 


A 


GRISM 


3 of 3 


15.01 


78.79 


0.68 


68.63 


used 


2007-04-14 


04:01:23.388 


P79 


C 


GRISM 


1 of 1 


43.49 


80.29 


0.68 


63.90 


used 


2007-04-14 


04:56:33.000 


P79 


B 


GRISM 


Oof 1 


24.46 


84.98 


0.68 


50.71 


not used 


2007-04-22 


23:36:20.000 


P79 


D 


GRISM 


lof8 


58.41 


58.54 


0.70 


63.41 


used 


2007-04-22 


03:48:12.779 


P79 


D 


GRISM 


2 of 8 


55.34 


81.89 


0.70 


58.22 


used 


2007-12-02 


07:58:13.000 


P80 


D 


PRISM 


Oof 3 


52.05 


48.04 


0.12 


50.30 


not used 


2007-12-02 


08:58:00.000 


P80 


D 


PRISM 


1 of 3 


58.65 


58.93 


0.12 


63.90 


used 


2007-12-12 


06:16:34.000 


P80 


F 


PRISM 


1 of 1 


71.54 


-16.21 


0.14 


36.30 


used 


2007-12-12 


07:28:52.000 


P80 


D 


PRISM 


2 of 3 


53.33 


50.29 


0.14 


52.60 


used 


2007-12-12 


08:19:04.000 


P80 


E 


PRISM 


1 of 1 


68.45 


126.5 


0.14 


64.00 


used 


2007-12-12 


09:05:20.166 


P80 


D 


PRISM 


2 of 3 


61.77 


64.26 


0.14 


74.50 


used 


2008-01-10 


07:13:33.000 


P80 


B 


PRISM 


2 of 2 


31.14 


65.29 


0.20 


76.50 


used 


2008-01-10 


05:03:32.000 


P79 


B 


GRISM 


1 of 1 


25.01 


44.22 


0.20 


45.49 


used 


2008-01-11 


04:54:29.000 


P80 


B 


PRISM 


1 of 1 


24.73 


43.11 


0.20 


45.60 


used 


2008-01-12 


05:27:07.000 


P80 


A 


PRISM 


Oof 2 


13.39 


50.69 


0.20 


52.60 


not used 


2008-01-12 


06:18:56.000 


P80 


A 


PRISM 


Oof 2 


14.65 


58.83 


0.20 


64.40 


not used 


2008-01-12 


04:37:00.000 


P79 


A 


GRISM 


Oof 1 


11.96 


39.65 


0.20 


41.26 


not used 


2008-02-20 


04:42:41.000 


P80 


D 


PRISM 


Oof 4 


62.75 


66.32 


0.28 


78.50 


not used 


2008-02-21 


03:27:22.000 


P80 


D 


PRISM 


Oof 3 


57.55 


57.17 


0.28 


61.20 


not used 


2008-02-22 


03:37:17.000 


P80 


B* 


PRISM 


Oof 6 


29.39 


59.09 


0.28 


64.30 


not used 


2008-02-22 


04:25:40.000 


P80 


B* 


PRISM 


Oof 6 


31.15 


65.28 


0.28 


75.30 


not used 


2008-03-11 


02:14:28.945 


P80 


A 


PRISM 


3 of 3 


14.29 


56.53 


0.31 


61.60 


used 


2008-03-11 


02:48:44.946 


P80 


A 


PRISM 


3 of 3 


15.05 


61.44 


0.31 


69.40 


used 


2008-03-13 


01:21:28.946 


P80 


B 


PRISM 


7 of 8 


26.17 


48.54 


0.32 


51.30 


used 


2008-03-13 


01:53:55.946 


P80 


B 


PRISM 


7 of 8 


27.93 


54.46 


0.32 


58.70 


used 


2008-03-13 


03:34:35.946 


P80 


B 


PRISM 


7 of 8 


31.64 


67.68 


0.32 


81.40 


used 


2008-03-25 


02:22:58.000 


P80 


C 


PRISM 


6 of 6 


46.74 


65.36 


0.34 


76.70 


used 


2008-03-25 


02:02:04.000 


P79 


C 


GRISM 


1 of 1 


45.74 


62.80 


0.34 


71.28 


used 


2008-03-27 


00:03:49.946 


P80 


A 


PRISM 


1 of 1 


12.56 


44.63 


0.34 


47.10 


used 


2008-03-27 


00:54:34.945 


P80 


C 


PRISM 


1 of 2 


41.87 


54.41 


0.34 


58.60 


used 


2008-04-02 


05:13:23.946 


P81 


C 


PRISM 


2 of 4 


40.95 


82.29 


0.35 


56.96 


used 


2008-04-07 


00:24:01.946 


P81 


A 


PRISM 


1 of 1 


14.18 


55.86 


0.36 


60.61 


used 


2008-04-07 


01:01:29.000 


P81 


C 


PRISM 


1 of 1 


45.24 


61.65 


0.36 


69.13 


used 


2008-04-07 


01:47:45.945 


P81 


B 


PRISM 


2 of 2 


31.47 


66.78 


0.36 


79.52 


used 


2008-04-08 


00:39:00.946 


P81 


A 


PRISM 


1 of 1 


14.64 


58.73 


0.36 


64.91 


used 


2008-04-08 


01:37:05.946 


P81 


B 


PRISM 


2 of 2 


31.32 


66.04 


0.36 


78.03 


used 


2008-04-28 


01:24:40.000 


P81 


D 


PRISM 


6 of 6 


63.91 


72.77 


0.40 


84.97 


used 


2008-11-30 


07:49:13.000 


P81 


D 


PRISM 


1 of 1 


51.02 


46.15 


0.81 


47.18 


used 


2008-12-27 


07:24:08.000 


P82 


A 


PRISM 


1 of 1 


14.84 


60.02 


0.86 


65.61 


used 


2008-12-27 


08:14:23.000 


P80 


C 


PRISM 


Oof 1 


47.04 


66.21 


0.86 


76.98 


not used 



...continues on next page 
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Date 


UT time 


pa 




Disp^ 


Car' 


(m) 


PA*' (°) 


Phase^ 




OF' 


2009-01-16 


05:01:58.946 


P82 


E 


PRISM 


Oof 1 


62.98 


124.5 


0.90 


57.84 


not used 


2009-01-16 


06:54:20.000 


P82 


F 


PRISM 


4 of 5 


71.35 


1.97 


0.90 


76.67 


used 


2009-01-19 


06:37:49.000 


P82 


D 


PRISM 


2 of 2 


62.41 


65.54 


0.91 


75.61 


used 


2009-01-19 


06:51:54.959 


P82 


D 


PRISM 


2 of 2 


62.80 


66.42 


0.91 


78.75 


used 


2009-01-20 


04:05:35.000 


P80 


C 


PRISM 


2 of 4 


36.00 


39.97 


0.91 


41.95 


used 


2009-01-20 


04:50:06.946 


P82 


C 


PRISM 


1 of 4 


39.53 


49.19 


0.91 


52.04 


used 


2009-01-20 


06:41:18.946 


P82 


B 


PRISM 


1 of 1 


31.23 


65.66 


0.91 


77.27 


used 


2009-01-22 


05:59:08.945 


P82 


B 


PRISM 


Oof 1 


30.12 


61.51 


0.91 


69.54 


not used 


2009-01-22 


06:07:04.959 


P82 


B 


PRISM 


Oof 1 


30.41 


62.53 


0.91 


71.34 


not used 


2009-01-22 


06:42:31.946 


P82 


C 


PRISM 


3 of 5 


47.18 


66.67 


0.91 


79.28 


used 


2009-01-25 


03:35:25.945 


P82 


A 


PRISM 


2 of 8 


11.64 


36.53 


0.92 


39.57 


used 


2009-01-25 


07:17:58.945 


P82 


A 


PRISM 


8 of 8 


16.01 


71.37 


0.92 


86.60 


used 


2009-02-19 


01:53:34.000 


P82 


D 


PRISM 


Oof 1 


46.81 


37.21 


0.97 


38.77 


not used 


2009-02-27 


06:54:56.945 


P82 


B 


PRISM 


1 of 1 


29.27 


79.89 


0.98 


64.50 


used 


2009-02-28 


03:28:46.946 


P82 


B 


PRISM 


2 of 2 


29.92 


60.86 


0.98 


68.43 


used 


2009-02-28 


04:04:29.946 


P82 


B 


PRISM 


2 of 2 


31.14 


65.27 


0.98 


76.49 


used 


2009-03-08 


00:42:31.015 


P82 


C 


PRISM 


3 of 3 


34.29 


34.38 


1.00 


37.82 


used 




09 -14-47 000 


P80 


Q 


PRISM 


3 of 3 


42.15 


55.01 


1.00 


58.75 


used 


2009-03-08 


02:52:36.000 


P80 


C 


PRISM 


3 of 3 


44.78 


60.63 


1.00 


67.36 


used 


2009-03-08 


03:50:00.000 


P82 


A 


PRISM 


1 of 1 


15.81 


67.46 


1.00 


80.24 


used 


2009-03-17 


03:16:00.000 


P82 


D 


PRISM 


2 of 2 


63.26 


67.61 


0.01 


80.54 


used 


2009-03-17 


03:48:12.000 


P82 


D 


PRISM 


2 of 2 


63.98 


70.78 


0.01 


86.38 


used 


2009-05-03 


23:43:46.000 


P83 


D 


PRISM 


Oof 2 


62.16 


65.02 


0.09 


75.35 


not used 



Notes. ^ Semester of proposal. ^ AT station: A = EO-GO, B = GO-HO, B* = AO-DO, C = EO-HO, D = DO-HO, E = DO-Gl or F = HO-Gl 
(cf. Table C). Dispersive element. ^ Calibrators used out of those available. ^ Projected baseline length. ^ Position angle of the projected baseline 
on the sky. § Visual light phase. ^ Elevation. ^ Quality flag showing if that observation is used for the model fitting or not. 
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